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Abstract. Soliton solutions of the KP equation have been studied since 1970, when Kadomtsev and 
Petviashvili proposed a two-dimensional nonlinear dispersive wave equation now known as the KP 
equation. It is well-known that one can use the Wronskian method to construct a soliton solution 
to the KP equation from each point of the real Grassmannian Gr^ n . More recently, several authors 
I3l ll4l l2ll4ll6l have studied the regular solutions that one obtains in this way: these come from points 
of the totally non-negative part of the Grassmannian (Gr^ n )>o- 

In this paper we exhibit a surprising connection between the theory of total positivity for the 
Grassmannian, and the structure of regular soliton solutions to the KP equation. By exploiting this 
connection, we obtain new insights into the structure of KP solitons, as well as new interpretations 
of the combinatorial objects indexing cells of (Gr^ , n )>o I23| . In particular, we completely classify 
the spatial patterns of the soliton solutions coming from (Gr^ n )>o when the absolute value of the 
time parameter is sufficiently large. We demonstrate an intriguing connection between soliton graphs 
for (Grj. n )>o and the cluster algebras of Fomin and Zelevinsky [9], and we use this connection to 
solve the inverse problem for generic KP solitons coming from (Gr^ , n )>o- Finally we construct all the 
soliton graphs for (Gr2, n )>o using the triangulations of an n-gon. 



The KP equation is a two-dimensional nonlinear dispersive wave equation which was proposed 
by Kadomtsev and Peviashvili in 1970 to study the stability problem of the soliton solutions of the 
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Korteweg-de Vries (KdV) equation [13] . The KP equation can also be used to describe shallow wa- 
ter waves, and in particular, the equation provides an excellent model for the resonant interaction of 
those waves (see [15] for recent progress). The equation has a rich mathematical structure, and is now 
considered to be the prototype of an integrable nonlinear dispersive wave equation with two spatial 
dimensions (see for example [22l [Tl 171 1211 [T2] ) . 

One of the main breakthroughs in the KP theory was given by Sato [25] , who realized that solutions 
of the KP equation could be written in terms of points on an infinite-dimensional Grassmannian. The 
present paper deals with a real, finite-dimensional version of the Sato theory; in particular, we are 
interested in solutions that are regular in the entire xy-plane, where they are localized along certain 
rays. Such a solution can be constructed from a point A of the real Grassmannian. More specifically, 
one can apply the Wronskian form [25, 26, 10, [12] to A to produce a T-function ta(x, y, t) which is a sum 
of exponentials, and from the T-function one can construct a solution ua(x, y, t) to the KP equation. 

Recently several authors have worked on classifying the regular soliton solutions [31 ® 01 [B] . These 
solutions come from points of the totally non-negative part of the Grassmannian, that is, those points 
of the real Grassmannian whose Plucker coordinates are all non-negative. They found a large variety of 
soliton solutions which were previously overlooked by those using the Hirota method of a perturbation 
expansion |12j . In the generic situation, the asymptotic pattern at y — > ±oo of the solution consists 
of n linc-solitons. However, because of the nonlinearity in the KP equation, the interaction pattern 
of the soliton solutions are very complex. Figure [T] illustrates the time evolution of the pattern of a 
soliton solution. Each figure shows the contour plot of the solution at a fixed time t in the xy-planc 
with x in the horizontal and y in the vertical directions. One of the main goals of this paper is to give 
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Figure 1. Time evolution of the spatial pattern of a soliton solution. The figures 



illustrate the contour plots of the solution (see Example 8.19 to reconstruct the figures) 



a combinatorial classification of such figures. 

Recently Postnikov [23J studied the totally non-negative part of the Grassmannian (Grk, n )>o from 
a combinatorial point of view. Total positivity has attracted a lot of interest in the last two decades, 
largely due to work of Lusztig [H][T9], who introduced the totally positive and non- negative parts of real 
reductive groups and flag varieties (of which the Grassmannian is an important example). Postnikov 
gave a decomposition of {Grk, n )>o into positroid cells, by specifying which Plucker coordinates are 
strictly positive and which are zero. He also introduced several remarkable families of combinatorial 
objects, including decorated permutations, d-diagrams, plabic graphs, and Grassmann necklaces, in order 
to index the cells and describe their properties. 

To understand a soliton solution ua(x, y,t) of the KP equation, one fixes the time t, and plots the 
contour plot C t (uA) in the xy-plane, the set of points where Ua(x, y, t) has a local maximum. This 
gives rise to a tropical curve in the plane; concretely, this shows the positions in the plane where the 
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corresponding wave has a peak, see Figure [l] Each region in the complement of the curve is naturally 
associated to a dominant exponential from the r-function TA{x,y,t) defined by ua — 2 ;|^2 In . We 
associate to Ct{uA) a soliton graph by forgetting the metric structure but remembering the dominant 
exponentials. 

In this paper we establish a tight connection between total positivity on the Grassmannian and the 
regular soliton solutions of the KP equation. This allows us to apply machinery from total positivity 
to understand soliton solutions of the KP equation. In particular: 

• we completely classify the soliton graphs coming from (G , r2 i „)>o, and classify the soliton graphs 
and contour plots coming from (Grk, n )>o when the absolute value of the time parameter t is 
sufficiently large, |i| > (Theorems [12. 1[ pD?] and |8. 12[ ) ; 

• we demonstrate an intriguing connection between soliton graphs for (Gr).,n)>o and the cluster 
algebras of Fomin and Zelevinsky [9] (Theorem 10.121; 

• we solve the inverse problem for KP solitons coming from (Grk, n )>o, and from (Grk, n )>o when 
\t\ » (Theorems [Il31|11.8| and |11.2[ ). 

In the other direction, we give a KP soliton interpretation to nearly all of Postnikov's combinatorial 



objects, as well as a new characterization of reduced plabic graphs (Theorem 10.51 



Note that prior to our work almost nothing was known about the classification of soliton graphs, 
except in the cases of (GV lin )>o [5], and (Gr2,4)>o [BJ. 

The structure of the paper is as follows. In Sections [2] and [3] we provide background on total 
positivity on the Grassmannian, and soliton solutions to the KP equation. In Section [4] we explain 
how to associate soliton graphs to soliton solutions of the KP equation. In the next four sections 
(Sections |5j [6j [7J and [8]) we explain the relationships between combinatorial objects labeling positroid 
cells and the corresponding soliton solutions. In particular, we explain how (decorated) permutations 
and Grassmann necklaces control the asympototics of soliton graphs when \y\ 3> 0, and how J-diagrams 
control the soliton graphs at \t\ ^> 0. We also explain the connection between plabic graphs and soliton 
graphs. In Section [9] we explain how the existence of A-crossings in contour plots corresponds to "two- 
term" Phicker relations. In Section [10] we prove that generically, the dominant exponentials labeling 
the regions of a soliton graph for (Grk, n )>o comprise a cluster for the cluster algebra of Gr^^ n . In 
Section 11 we address the inverse problem for regular soliton solutions to the KP equation. Finally, in 
Section 12 we completely classify the soliton graphs coming from (Gr2 jn )>o, and construct them all 



using triangulations of an n-gon. 

The present paper provides proofs of the results announced in |16j . In the sequel to this work |17j we 
have extended many of the results of the present paper from the non-negative part of the Grassmannian 
to the real Grassmannian. In a future paper we plan to make a detailed study of the relationship between 
cluster transformations and the evolution of soliton graphs. 

Acknowledgements: The first and second authors are grateful for the hospitality of the math 
departments at UC Berkeley and Ohio State, where some of this work was carried out. They are also 
grateful to Sara Billey, whose comments helped them to greatly improve the exposition. 



2. Total positivity for the Grassmannian 

In this section we review the Grassmannian Grk, n and Postnikov's decomposition of its non-negative 
part (Grfc jn )>o into positroid cells [23]. Note that our conventions slightly differ from those of [23] . 

The real Grassmannian Gr^ n is the space of all fc-dimensional subspaces of E™. An element of Gr^ n 
can be viewed as a full-rank k x n matrix modulo left multiplication by nonsingular k x k matrices. In 
other words, two k x n matrices represent the same point in Gr^.n if and only if they can be obtained 
from each other by row operations. Let ('^) be the set of all fc-element subsets of [n] :={!,..., n}. For 
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/ G ( t ), let A i (A) denote the maximal minor of a k x n matrix A located in the column set /. The 
map A i — y (Ai(A)), where / ranges over (^'), induces the Pliicker embedding Gr^.n KPw -1 . 

For Ai C ( the matroid stratum SVi is the set of elements of Gr^ n represented by all k x n 
matrices A with A/ (A) ^ for / e and Aj(A) — for J ^ Ai. The decomposition of Gr^ n into 
the strata SVi is called the the matroid stratification. 

Definition 2.1. The totally non-negative Grassmannian (Gri~ t n)>o (respectively, totally positive Grass- 
mannian (Grk t n)>a) is the subset of Gr^^ that can be represented by k x n matrices A with all Ai(A) 
non- negative (respectively, positive). 

Postnikov [22] studied the decomposition of (GVfc )n )>o induced by the matroid stratification. More 
specifically, for Ai C (^), he defined the positroid cell SjJJ™ as the set of elements of (Grk. n )>o 
represented by all k x n matrices A with Ai(A) > for / e and Aj(A) — ioi J ^ Ai. It turns 
out that each nonempty S 1 ^™ is actually a cell [23], and that this decomposition of (Grk, n )>o is a CW 
complex [24]. Note that (Grk, n )>o is a positroid cell; it is the unique positroid cell in (Grk, n )>a 01 
top dimension k(n — k). Postnikov showed that the cells of (Gr^ , n )>o are naturally labeled by (and in 
bijection with) the following combinatorial objects [23] : 

• Grassmann necklaces I of type (fc, n) 

• decorated permutations tt on n letters with k weak excedances 

• equivalence classes of reduced plabic graphs of type (k, n) 

• J-diagrams of type (k,n). 

For the purpose of studying solitons, we are interested only in the irreducible positroid cells. 

Definition 2.2. We say that a positroid cell is irreducible if the reduced- row echelon matrix A of 
any point in the cell has the following properties: 

(1) Each column of A contains at least one nonzero element. 

(2) Each row of A contains at least one nonzero element in addition to the pivot. 

The irreducible positroid cells are indexed by: 

• irreducible Grassmann necklaces 1 of type (fc, n) 

• derangements jronn letters with k excedances 

• equivalence classes of irreducible reduced plabic graphs of type (fc, n) 

• irreducible J-diagrams of type (fc,n). 

We now review the definitions of these objects and some bijections among them. 

Definition 2.3. An irreducible Grassmann necklace of type (fc, n) is a sequence X = . . . , I n ) of subsets 
I r of [n] of size k such that, for i £ [n], = (L L \ {i}) U {j} for some j ^ i. (Here indices i are taken 
modulo n.) 

Example 2.4. 1 = (1257,2357,3457,4567,5678,6789,1789,1289,1259) is an example of a Grassmann 
necklace of type (4,9). 

Definition 2.5. A derangement tt — (tv(1), . . . , 7r(n)) is a permutation n £ S n which has no fixed points. 
An excedance of 7r is a pair (i,7r(i)) such that > i. We call i the excedance position and ir(i) the 
excedance value. Similarly, a nonexcedance is a pair (i, 7r(i)) such that ir(i) < i. 

Remark 2.6. A decorated permutation is a permutation in which fixed points are colored with one of two 
colors. Under the bijection between positroid cells and decorated permutations, the irreducible positroid 
cells correspond to derangements, i.e. those decorated permutations which have no fixed points. 

Example 2.7. The derangement tt = (6, 7, 1, 2, 8, 3, 9, 4, 5) G Sg has excedances in positions 1, 2, 5, 7. 



KP SOLITONS AND TOTAL POSITIVITY FOR THE GRASSMANNIAN 



5 



Definition 2.8. A plabic graph is a planar undirected graph G drawn inside a disk with n boundary 
vertices 1, ...,n placed in counterclockwise order around the boundary of the disk, such that each 
boundary vertex i is incident to a single edgeQ Each internal vertex is colored black or white. See 
Figure [2] for an example. 
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k = 6, n = 16 

^ = (10,9,9,8,5,2) 



Figure 2. A plabic graph, and an irreducible Le-diagram L = (A,D)/ C) 



Definition 2.9. Fix k and n. Let Y\ denote the Young diagram of the partition A. A J-diagram (or 
Le-diagram) L = (A, D)^ ^ of type (fc,n) is a Young diagram Y\ contained in a k x (n — k) rectangle 
together with a filling D : Y\ — > {0, +} which has the d-property: there is no which has a + above it 
in the same column and a + to its left in the same row. A J-diagram is irreducible if each row and each 
column contains at least one +. See the right of Figure [2] for an example of an irreducible J-diagram. 

Theorem 2.10. [23 Theorem 17.2] Let Sffi be a positroid cell in (Gr k ,„)>o. For 1 < r < n, let I r be 
the element of M which is lexicographically minimal with respect to the order r<7' + l<---<n< 
1 < 2 < . . . r — 1. Then := (Ii, . . . , /„) is a Grassmann necklace of type (k, n). 

Lemma 2.11. [23l Lemma 16.2] Given an irreducible Grassmann necklace I, define a derangement 
7r = 7r(Z) by requiring that: if — (Ii \ {i}) U {j} for j 7^ i, then n(j) = i^ Indices are taken modulo 
n. Then I — > 7r(I) is a bijection from irreducible Grassmann necklaces I = (/]_,..., /„) of type (fc, n) 
to derangements tt(Z) S S n with k excedances. The excedances of 7r(I) are in positions 1%. 



Example 2.12. If X and 7r are defined as in Examples 2.4 and 2.7 then 7r(I) 



Definition 2.13. Given a J-diagram L contained in a k x (n — k) rectangle, label its southeast border 
with the numbers 1,2, ... ,n, starting at the northeast corner. Replace each + with an "elbow" and 
each with a "cross"; see Figure [3] Now travel along each "pipe" from southeast to northwest, and 
label the end of a pipe with the same number that labeled its origin. Finally, we define a permutation 
7r = tt(L) as follows. If i is the label of a vertical edge on the southeast border of L, then set ir(i) equal 
to the label of the vertical edge on the other side of that row. If i is the label of a horizontal edge on 
the southeast border of L, then set equal to the label of the horizontal edge on the opposite side 
of that column. See Figure [3j 



Proposition 2.14. The map defined above gives a bijection from irreducible J-diagrams contained in a 
k x (n — k) rectangle to derangements on n letters with k excedances. 

Proof. This map can be shown to coincide with that from |28l Section 2], and up to a convention 
change, coincides with the map in [2"31 Corollary 20.1]. □ 



The convention of 1231 was to place the boundary vertices in clockwise order. 

'Postnikov's convention was to set 7r(i) = j above, so the permutation we are associating is the inverse one to his. 
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FIGURE 3. The Le-diagram L together with the computation of ir(L) = (7, 4, 2, 9, 1, 3, 8, 6, 5). 

Remark 2.15. Consider a positroid cell Sffi, and suppose that the Grassmann necklace 1, the derange- 
ment 7r, and the J-diagram L, satisfy 1 = 1(A4), and ir — tt(1) — tt(L). Then we also refer to this cell 

00 Qtnn Qtnn Qtnn Qtnn Qtnn 

db °M(I)' D M(L)' °I ' J ir ' D L j clt. 

3. SOLITON SOLUTIONS TO THE KP EQUATION 

We now explain how to obtain a soliton solution to the KP equation from a point of {Gr^ n )> . 

3.1. From a point of the Grassmannian to a r- function. We start by fixing real parameters k.- l 

v 

such that K\ < Ki ■ ■ ■ < k„, which are generic, in the sense that the sums Ki. are all distinct for 

i=i 

any p with 1 < p < n. We also assume that the differences between consecutive k^s are similar, that 
is, Ki+i — Ki is of order one. 

Let {Ei \ i — 1, . . . , m} be a set of exponential functions in (ti, . . . ,t m ) E M. n defined by 

Ei(ti, ...,t m ):= exp i^tA . 

Here m is an arbitrary positive integer which is at least 3. Thinking of each Ei as a function of t\, we 
see that the elements {Ei} are linearly independent, because their Wronskian determinant with respect 
to ti is non zero: 

Wr(E u ...,£„):= det[(^'- 1) ) 1 < JJ <„] = - m) E x ---E n ± 0. 

Here E\ j) := d 3 E t /dt[ = n\E t . 

Let A be a full rank k x n matrix. We define a set of functions {/i, . . . , /&} by 

(/l>/2; ■ ■ ■ ,fk) T = A ■ (E 1: E 2 , ■ ■ ■ ,E n ) T , 

where (. . .) T denotes the transpose of the vector (. . .). 

Since the exponential functions {Ei} are linearly independent, we identify them as a basis of R™, and 
then {/i, . . . , fk} spans a fc-dimensional subspace. This identification can be seen, more precisely, as 
Ei <H- (1, Ki, . . . , k™ _1 ) t € IR". This subspace depends only on which point of the Grassmannian Grk. n 
the matrix A represents, so we can identify the space of subspaces {f%, . . . , /&} with Grk, n - 

The T-function of A is defined by 

(3.1) T A (h,...,t m ) :=W(fi,f 2 ,...J k ). 

For / = {ii,...,i k } e we set 

Ei(tx,...,t m ) :=Wr(E il ,E i2 ,...,E ik ) = JJ (n im - m e ) E i± ■ ■ ■ E ik > 0. 
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Applying the Binet-Cauchy identity to the fact that fj = ^ CbjiEi, we get 
(3.2) T A (t u ...,t m )= J2 A I (A)E I (t 1 ,...,t m ). 

It follows that if A £ (Gr k , n )>0, then t a > for all (ti, . . . , t m ) € R m . 



Thinking of t a as a function of A, we note from (3.2) that the r-function encodes the information 
of the Pliicker embedding. More specifically, if we identify each function Ej with I = . . . , i^} with 
the wedge product E il A • • ■ A E ik (recall the identification E^ o (1. Hi, ... , k™ _1 ) t ), then the map 

4 ^ ta has the Pliicker coordinates as coefficients. 
3.2. Prom the r-function to solutions of the KP equation. The KP equation for u(x,y,t) 

d ( ^du du d 3 u\ 2,® 2u (J 
dx \ dt dx dx 3 J dy 2 

was proposed by Kadomtsev and Petviashvili in 1970 [13] . in order to study the stability of the soliton 
solutions of the Korteweg-de Vries (KdV) equation under the influence of weak transverse perturbations. 
The KP equation can be also used to describe two-dimensional shallow water wave phenomena (see for 
example |15j). This equation is now considered to be a prototype of an intcgrablc nonlinear partial 
differential equation. For more background, see [331 E E21 HI] • 

It is well known (see [12l 01 El [6]) that if we set x = t\, y = t 2 , t — t 3 (treating the other t^s as 



constants), then the r-function defined in (3.1 1 provides a soliton solution of the KP equation, 



d 2 

(3.3) u A {x,y,t) = In ta (x, y,t). 

The parameters tj correspond to the symmetries of the KP equation, and the set of those flows defined 
by the symmetries is called the KP hierarchy. 

It is easy to show that if A £ (Grk. n )>a, then such a solution ua(x, y, t, t±, . . . , t m ) is regular for all 
tj e K. In the sequel to this paper [T7], we will show that if UA(x,y,t) is regular for all x,y,t, then 
A e (GVfc.„)>o- (In [H Proposition 4.1], a weaker statement was proved: if UA(x,y,t,t4, . . . ,t m ) is 
regular for all tj € K, then A 6 (GVfc,n)>o-) F° r this reason we are interested in solutions UA{x,y,t) of 
the KP equation which come from points A of (GVfc n )>o. Throughout this paper when we speak of a 



soliton solution to the KP equation, we will mean a solution UA(x,y,t) which has form (3.3) 



4. From soliton solutions to soliton graphs 

Here we explain the fundamental structure of the soliton solution and define the contour plot of the 
soliton solution and the corresponding soliton graphs. 

4.1. Contour plots and soliton graphs. One can visualize a solution ua{x, y,t) in the xy-plane by 
drawing level sets of the solution when the coordinates (t, <4, . . . , i m ) are fixed. For each r e I, we 
denote the corresponding level set by 

C r (t, t 4 ,... 1 t m ) := {{x, y) e M 2 : u A (x, y, t,U,..., t m ) = r}. 

Figure [4] depicts both a three-dimensional image of a solution ua{x, y, t), as well as multiple level sets 
C r . These levels sets are lines parallel to the line of the wave peak. 

Example 4.1. We compute the soliton solution Ua(x, y,t) associated to the matrix A = (1 a) with 
a > 0, considered as an element of (Gr l!2 )>o- Write E\ = e Bl and aE 2 = e® 2 . Then the r-function ta 
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and the soliton solution ua are given by 

T A (x,y,t) = e 01 +e 92 = 2e5 (ei+02) coshi(6»i - 2 ), and 

u A (x,y,t) = - K 2 ) 2 secri 2 i(6>i - 2 ). 

This is a line-soliton solution, which represents a two-dimensional wave. Since sechx takes its maximum 
value at x = 0, the wave crest is on the line 0% — 02, i-e. 

x + (kj + n-))y + {k\ + K\K2 + K 2 )^ + constant = In a. 

K 2 - Kl 




Figure 4. A line-soliton solution ua(^, y,t) where A = (1,1) G (Gri ]2 )>o, depicted 
via the 3-dimensional profile ua(x, y, 0), and the level sets of ua{x, y, 0). E$ represents 
the dominant exponential in each region. 



If we fix t, the line Q\ = 62 divides the xy-plane into two regions. The exponential E\ dominates 
in the region including x <C 0, and E 2 dominates the other region where i>0, We label each region 
by its dominant exponential. Figure [4] depicts Ua(x, y,t), where t = t± = ■ ■ ■ = t m = 0, a = 1, and 
(«i,«a) = (-l,2). 

To study the behavior of ua(x, y, t) for A G S/Jj" an< ^ fi xm § all other time variables t^, . . . , t m , we 
consider the dominant exponentials at each point (x,y,t), and we define 

f A (x,y,t) = max{Aj(A)Ej(x,y,t,t4,, . . . ,t m )} = max{Aj(A)KjE n ...E 3k } 
= max|exp^ln(A J (A)X J ) + f^K^x + n\y + K %t + o) J. 



where Kj :— Yli<m( K 3m ~ K h) > and Q '•••> m — J2 From (3.2), we see that ta can be 



approximated by /a- 

Let fA{x,y,t) be the closely related function 



(4.1) fA{x,y,t) = max {ln(Aj(A)Kj) + £ (it jt x + nly + «f.t + i 1 '-' m )\ . 

JGM [ j_j J * Jl Jz J 

Note that at a given point (x, y, t), Ja{x, y, t) is equal to a given term if and only if f A (x, y, t) is equal 
to the exponentiated version of that term. 
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Definition 4.1. Given a solution UA{x,y,f) of the KP equation as in (3.3), we define its contour plot 
C{ua) to be the locus in R 3 where fA(x,y, t) is not linear. If we fix t = t , then we let C to (uA) be the 
locus in R 2 where Ja(x, y,t = to) is not linear, and we also refer to this as a contour plot. 

Remark 4.2. In the majority of this paper, we will assume that ti = ts = • • • = t m = 0. However, in 



Section 12 we will use these extra variables in a non-trivial way. 



Remark 4.3. C ta { u A) provides an approximation of the location of the wave crests. Here the time to 
should be generic in the sense that the number of regions in the plot Ct (A) remains the same as those 
in Ct(uA) for the time interval \t — to\ = 0(1). 



It follows from Definition 4.1 that C(ua) and Ct (v,A) are piecewise linear subsets of IR and R , 
respectively, of codimension 1 . 

Proposition 4.4. If each Kj is an integer, then C(ua) is a tropical hypersurface in ]R 3 , and C to {uA) is a 
tropical hypersurface (i.e. a tropical curve) in R 2 . 



Proof. This follow immediately from Definition 4.1 See [TTJ Section 1.1] or [2^ for an introduction 
to tropical hypersurfaces. □ 

Note that each region of the complement of C(ua) and C to (ua) hi K 3 and R 2 , respectively, is a domain 
of linearity for fA(x,y,t), and hence each region is naturally associated to a dominant exponential 



Aj(A)Ej(x,y,t) from the r-function (3.2 1. We label this region by J or Ej. 

The contour plot Ct {uA) consists of line segments, some of which have finite length, while others 
are unbounded and extend in the y direction to ±oo. We call these unbounded lines asymptotic or 
unbounded linc-solitons. Among the finite line segments, some of them represent phase shifts and have 
lengths which are determined by the K-parameters (see \§\ for details). We will ignore the phase shifts in 
this paper: they are negligable, provided that each Kj+i — Ki is of order 1, and that we consider contour 
plots with large scales of (x,y). We call the other finite line segments simply line-solitons. Note that 
each line-soliton represents a balance between two dominant exponentials in the r-function. 

The following lemma is based on the assumption that each line separating two regions in the contour 
plot is a line-soliton, that is, we ignore the phase shifts. 

Lemma 4.5. [SJ Proposition 5] The index sets of the dominant exponentials of the r-function in adjacent 
regions of the contour plot in the xy-plane are of the form {i, TO2, . . . , m^} and {j, rri2, ■ ■ . , rrik}. 



According to Lemma |4.5[ those two exponential terms have k — 1 common phases, so we call the line 
separating them a line-soliton of type Locally we have 



t a « A I (A)E I + Aj(A)Ej = (A I (A)K I E i + Aj(A)K jEj) J[ E m „ 

1=2 

Aj(A)Kj 



so the equation for this line-soliton is 

(4.2) x + (Ki + Kj)y + (k 2 + KiKj + n 2 At + constant — In 



A I (A)K I 



Note that the slope of the line is determined by Ki + Kj, and the constant on the left-hand side of (4.2) 
is if t& = ■ ■ ■ = t m = 0. In that case the location of the line is determined by the ratio of the Pliicker 
coordinates corresponding to the dominant exponentials on either side of the line-soliton. 



Remark 4.6. Consider a line-soliton given by (4.2). Compute the angle ^uj] between the line-soliton 
of type [i,j] and the positive y-axis, measured m the counterclockwise direction, so that the negative 
a;-axis has an angle of w and the positive a;-axis has an angle of — ^. Then tan^^ j] = k 
refer to K t + Kj as the slope of the [i, j] line-soliton (see Figure [2]). 
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We will be interested in the combinatorial structure of a contour plot, that is, the pattern of how 
line-solitons interact with each other. 

Generically we expect a point at which several line-solitons meet to have degree 3; we regard such 
a point as a trivalent vertex. Three line-solitons meeting at a trivalent vertex exhibit a resonant 



interaction (this corresponds to the balancing condition for a tropical curve). See Section 4.2 One may 
also have two line-solitons which cross over each other, forming an X-shape: we call this an X-crossing, 
but do not regard it as a vertex. See Figure[5]for examples. As mentioned above, in general, there exists 
a phase-shift at each X-crossing (see [B]). However the line-segments corresponding to the phase-shifts 
have the constant lengths, and we ignore them in this paper. Vertices of degree greater than 4 are also 
possible. 

Definition 4.7. A contour plot C t (uA) is called generic if all interactions of line-solitons are at trivalent 
vertices or are X-crossings. 

The following notion of soliton graph forgets the metric data of the contour plot, but preserves the 
data of how line-solitons interact and which exponentials dominate. 

Definition 4.8. Let Ct {uA) be a generic contour plot with n unbounded line-solitons. Color a trivalent 
vertex black (respectively, white) if it has a unique edge extending downwards (respectively, upwards) 
from it. Label a region by Ej if the dominant exponential in that region is AjEj. Label each edge 
(line-soliton) by the type of that line-soliton. Preserve the topology of the metric graph, but forget 
the metric structure. Embed the resulting graph with bicolored vertices into a disk with n boundary 
vertices, replacing each unbounded line-soliton with an edge that ends at a boundary vertex. We call 
this labeled graph the (generic) soliton graph Gt a {uA)- 

Abusing notation, we will often refer to the edges of Gt {uA) as line-solitons, and use the terminology 
unbounded line-solitons and unbounded regions to refer to the edges and regions incident to the boundary 
of the disk. 




Figure 5. Example of a contour plot and soliton graph associated to S^ 1 " 1 with 
7r= (7,4,2,9,1,3,8,6,5). 

See Figure [5] for an example of a soliton graph. Although we have not labeled all regions or all edges, 
the remaining labels can be determined using Lemma |4.5| 

4.2. Resonance of solitons. Here we explain the physical meaning of the trivalent vertex in the 
contour plot. Recall that a line-soliton of [i, j]-type has the form 

u = ~{ni - Kjfsech^K^-] • x + fi^t), 
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where Kuji and are called the wavenumber-vector and the frequency, and they are given by 

(4.3) K [id] := {Kf l]V K y {i]] ) = ( Kj - m, k) - «?), %,] = Kj 3 - 4 . 

One should note here that there is an algebraic relation, called the dispersion relation of the KP equation, 
among Kuji and Quj], which is given by 

(4-4) D(K [i:j] ,Q [id] ) := -An^Kfa + (K^f + 3(i^,) 2 = 

This implies that if a plane wave of the form <fi(Kj -x+fljt) with the wavenumber-vector Kj = (KJ, K y ) 
and the frequency flj is a solution of the KP equation, then Kj and f2j must satisfy the dispersion 



relation. Note that the wavenumber- vectors and the frequency given in the form (4.3) are the solution 



of the algebraic equation (4.4), i.e. _D(Kj,f2j) = 0. 

In the wave theory in general, if the dispersion relation D(K, £1) = has the following property with 
two plane waves, 0i(Kj ■ x + flit) for i = 1,2, 

D(K 1 +K 2 ,n 1 +Q 2 ) = 0, 

then as a result, a third wave can be generated and the new wave ipsCK^,^) satisfies the so-called 
resonant conditions, 

K 3 = Ki+K 2 , and fi 3 = fii+0 2 . 

In the KP dispersion relation, the line-solitons of of types [j, P\, and [i,£] (here i < j < £) trivially 

satisfy the resonant conditions, i.e. 

(4.5) K[ ir £j = K^j] + Ky^j, = %j] + &\j,t]- 

This means that the resonant interactions happen quite naturally in the KP equation, and each 3-wave 
resonant interaction appears as a trivalent vertex in the contour plot. Then at the trivalent vertex, since 
the slope of each soliton is given by tan 'I'u-ji = Ki + Kj, those three solitons appear as [i,j], and 
[i, I] in counterclockwise order. This condition provides a key characterization of reduced plabic graphs 



in Section 10 Note that equation (4.5) corresponds to the balancing condition for a tropical curve. 



5. Permutations and soliton asymptotics 

Given a soliton graph Gt (uA) where A belongs to an irreducible positroid cell, we show that the 
labels of the unbounded solitons allow us to determine which positroid cell A belongs to. Conversely, 
given A in the irreducible positroid cell S*™ , we can predict the asymptotic behavior of the unbounded 
solitons in C to (Uj\). 

Theorem 5.1. Suppose A is an element of an irreducible positroid cell in (Grk, n )>o- Consider the 
contour graph Ct a (uA) for any time t$. Then there are k unbounded line-solitons at y 3> which are 
labeled by pairs [e r ,j r ] with e r < j r , and there are n — k unbounded line-solitons at y <C which are 
labeled by pairs [« r ,g r ] with i r < g r . We obtain a derangement in S n with k excedances by setting 
7r(e r ) = j r and Tr(g r ) = i r - Moreover, A must belong to the cell S'£ n ™. 

The first part of this theorem follows from work of Biondini and Chakravarty [2, Lemma 3.4 and 



Theorem 3.6] (see Proposition 5.2 below) and Chakravarty and Kodama [H Prop. 2.6 and 2.9], [SJ 
Theorem 5] (see Theorem 5.3 below). In particular, Chakravarty and Kodama had already associated 
a derangement n to A, but it was not clear how this it was related to the derangement indexing the cell 
containing A. Our contribution is a proof that th e de rangement ir is precisely the derangement labeling 
the cell S% m that A belongs to (see Proposition |5.4[ p| This fact is the first step towards establishing 
that various other combinatorial objects in bijcction with positroid cells (Grassmann necklaces, plabic 
graphs) carry useful information about the corresponding soliton solutions. 

3 S. Chakravarty informed us that he also proved an equivalent proposition. 
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Given a matrix A with n columns, let A(k, ...,£) be the submatrix of A obtained from columns 
k, + . . . , I— l,t, where the columns are listed in the circular order k, fc + 1, . . . ,n — 1, n, 1, 2, . . . , k— 1. 

Proposition 5.2. Lemma 3.4] Let A be a k x n matrix representing an element in an irreducible 
positroid cell in (GVfc ]ra )>o, and consider the contour plot Ct (uA) for any time t. Then there are n — k 
unbounded line-solitons at y <C and k unbounded line-solitons at y 0: 

There is an unbounded line-soliton of C t (u^) at y <C labeled [i,g] with i < g if and only if 

(5.1) rankv4(i, ... ,17 — 1) = rankA(z + 1,. .. ,g) — r&nkA(i, . . . ,g) — rank^4(i + 1, . . . , g — 1) + 1. 

Moreoever, g is a non-pivot column of A. 

And there is an unbounded line-soliton of Ct(uA) at y S> labeled [e, j] with e < j if and only if 

(5.2) rank A(j, . . . , e - 1) = rankA(j + 1, . . . , e) = rank A(j, . . . , e) = rank^(j + 1, . . . , e - 1) + 1. 
Moreoever, e is a pivot column of A. 

Theorem 5.3. [4, Prop. 2.6 and 2.9] [6 S Theorem 5] Consider an irreducible positroid cell Sjff 1 in 
(Grfc in )>o, and let A be a full rank matrix representing a point in that cell. Use the notation of 



Proposition 5.2 Define it := tt-(A4) by setting 7r(e) :— j and n'(g) := (i) for each pivot e and 
non-pivot c 

Proposition 

7T ! (M) =7T 



non-pivot 5. Then 7r is a derangement on n letters with k weak excedances 

7 A4 



Proposition 5.4. Consider an irreducible positroid cell S^ n = S 1 ^"™, where 7r = ir(I(M)). Then 



Proof. Consider a k x n matrix A representing an clement in S/\ n . Then all maximal minors of A are 
non- negative, and the column indices of the non-zero minors are the subsets in M.. Let us first consider 
the derangement 7r = tt(T(A4)). Let Ii = {i = x\,xi, ■ ■ ■ ,%k} be the lexicographically minimal minor 
in M. with respect to the total order i < i + 1 < ■ ■ • < n < 1 < ■ ■ ■ < i — 1. Then = (Ii \ {i}) U {j} 
is obtained from Ii\{i} by considering the column indices in the order i+l,i + 2,...,n, 1,2, ... ,i and 
greedily choosing the earliest index h such that the columns of A indexed by the set {X2, . . . , Xfe} U {h} 
are linearly independent. Then 7r(/i) is defined to be i. 

Now consider the ranks of various submatrices of A obtained by selecting certain columns. 

Claim 0. rank ^4(i + 1, . . . , h — 1, h) = 1 + rank A(i + 1, . . . , h — 1). This claim follows from the way 
in which we chose h above. 

Claim 1. rank y4(i, . . . , h) = rank A(i, . . . , h— 1). To prove this claim, we consider two cases. 
Either x\ <; h <, Xk or x\ <i Xk <i h, where <i is the total order i < i + 1 < • • • < n < 1 < • • • < i — 1. 
In the first case, the claim follows, because h is not contained in the set Ii but is contained in /j+i- In 
the second case, rank A(i, i + 1, i + 2, . . . , x^) = fc, and the index set {i, i + 1, . . . , 2^} is a strict subset 
of + 1, ... , /i}, so rank^4(i, . . . , h) = rank A(i, . . . , h — 1) = k. 

Now let R = rank J 4(z + 1, i + 2, . . . , h - 1). By Claim 0, rankA(i + 1, . . . , h) = R + 1. Therefore we 
have rankA(i, . . . , h) > rank A(i+1, . . . , h) = R+l. By Claim 1, rankA(i, . . . , h) = ra.nkA(i, . . . , h— 1), 
but rank A(«, . . . , h — 1) < ii + 1, so rank ^4(i, . . . , h) < R + 1. We now have rank A(i, . . . , h) = R+l. 
But also rank^4(i, . . . , h — 1) = ra,nkA(i, . . . , h) = R + 1. 

We have just shown that rank A(i, i + 1, . . . , h — 1) = rank A(i + 1, . . . , h — 1, h) = rank A(i, . . . , h) = 
rankA(i + 1, . . . , h — 1) + 1. Comparing these rank conditions to either part of Proposition |5.2| and 



using Theorem 5.3 we see that Tt'(h) — i. This shows that tt and 7r coincide. □ 

We now give a concrete algorithm for writing down the asymptotics of the soliton solutions of the 
KP equation. 

Theorem 5.5. Fix real generic parameters Ki < • • • < n n . Let A be a point in an irreducible positroid 
cell S*™ 71 in (Gr^„)> . (So n has k excedances.) For any to, the asymptotic behavior of the contour 
plot C to (uA) - its unbounded line-solitons and the dominant exponentials in its unbounded regions - 
can be read off from 7r as follows. 
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FIGURE 6. Asymptotic line-solitons for ir = (6,7,1,2,8,3,9,4,5). Each Eijki shows 
the dominant exponential in this region. 



For !/>0, there is an unbounded line-soliton which we label [i, for each excedance 7r(i) > i. 
From left to right, list these solitons in decreasing order of the quantity 

For 1/ < 0, there is an unbounded line-soliton which we label for each nonexcedance 

7r(j) < j. From left to right, list these solitons in increasing order Kj + k^n)- 
Label the unbounded region for i«0 with the exponential -E^,....^ , where i\, . . . ,ik are the 
excedance positions of tt. 



Use Lemma 4.5 to label the remaining unbounded regions of the contour plot. 



Proof. The fact that the set of unbounded line-solitons are specified by the derangement n comes 



from [SJ Theorem 5, page 125] and [51 Corollary 1, page 124]. It then follows from Remark 4.6 that for 
sufficiently large y (respectively, sufficiently small y), these solitons are ordered from left to right by 
decreasing (respectively, increasing) order of their slopes Ki + Kj. □ 

Example 5.6. Consider the positroid cell corresponding to (6, 7, 1, 2, 8, 3, 9, 4, 5) <E Sg. The algorithm of 
Theorem 5J3 gives rise to the picture in Figure [6j Note that if one reads the dominant exponentials in 



counterclockwise order, starting from the region at the left, then one recovers exactly the Grassmann 
necklace from Examples |2.4| and |2.12| This correspondence will be generalized in Proposition |5.4| and 
Theorem 16.21 



6. Grassmann necklaces and soliton asymptotics 

One particularly nice class of positroid cells is the TP or totally positive Schubert cells. A TP Schubert 
cell is the intersection of a usual Schubert cell with (Gr^. „)> . When A lies in a TP Schubert cell S^™", 
we can make another link between the soliton solution ua(%, y,t) and the combinatorics of (Grk, n )>o- 
Namely, the dominant exponentials labeling the unbounded regions of the soliton graph Gt(uA) form 
the Grassmann necklace associated to S*™ . 

It is easy to verify the following lemma. 

Lemma 6.1. A positroid cell S t J in = S'f nrl of {Grk, n )>o is a TP Schubert cell if and only if either 
(equivalently, both) of the following conditions hold: 

• The J-diagram L has each box filled with a +. 

• If i\ < i 2 < ■ ■ ■ < ik and j\ < 32 < •■ < jn-k are the positions of the excedances and 
nonexcedances, respectively, of tt, then 7r(«i) = n — k + 1, 7r(i2) = n — k + 2, . . . , 7r(ife) = n and 
tO'i) = 1 > n (h) = 2, ■ • ■ , Tr(jn-k) =n-k. 



We have the following result. 
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[ik,Jt(ik)] [i2,Jt(i2)] [h,jt(h)] 




Figure 7. Grassmann necklace R2, . . . , R n ) in a soliton graph associated to the 
TP Schubert cell S^ n . 



Theorem 6.2. Let A be an element of a TP Schubert cell Sl nn , and consider the soliton graph Gt {uA) 
for an arbitrary time to. Let the index sets of the dominant exponentials of the unbounded regions 
of G to {uA) be denoted i?i,...,i?„, where R\ labels the region at x <C 0, and i? 2 ,...,i?„ label the 
regions in the counterclockwise direction from J?x- Then (R\, . . . , R n ) is a Grassmann necklace X, and 

7r(Z) = 7T. 



Theorem 6.2 is illustrated in Example 5.6 See also Figure [6} 



Remark 6.3. Theorem 6.2 does not hold if we replace "TP Schubert cell" by "positroid cell." For exam- 
ple, the Grassmann necklace associated to the derangement it = (4, 3, 1, 2) is (12, 23, 34, 24). However, if 
Ki — 0, «2 = 1, K3 — 1.5, K4 = 1.75, then the corresponding sequence of dominant exponentials labeling 
the unbounded regions of any soliton graph coming from the cell is (12, 23, 34, 13). 



Remark 6.4. To recover a Grassmann necklace I = (Ji, 
the procedure of Lemma 2.11 ), we do the following: 



, I n ) from a derangement ir G S n (inverting 



• Set Ii = {ii, . . . , ifc}, the positions of the excedances of n. 

• For each r > 1, set I r+1 = (I r \ {r}) U {^(r)}. 



We now prove Theorem 6.2 
Proof. 



Let ix < 



be the positions of the excedances of n, and let ji < ■ ■ ■ < j n -k be 
the positions of the nonexcedances. By Lemma 6.1 we have that 7r(ii) < nfa) < ■ ■ ■ < 7r(ifc), and 
7r (j 1 ) < 7r(j 2 ) < • • • < 7r(jn-fc)- Define the partial order -< on pairs of integers in {1,2,..., n} by 
setting -< if and only if Ki + Kj < Ki> + Ky. Then the condition that Kx < K2 < ■ ■ ■ < n n 

implies that (ii,7r(ii)) -< (« 2 ,7r(z 2 )) -<•••-< (ik,^(ik)) and (ji,7r(ji)) -<•••-< (j n -k,^(Jn-k))- Using 



Theorem 5.5 the asymptotic directions of the contour graph of the soliton solution are as in Figure [7| 
From the conditions on our permutation, we must have tt(ji) = 1, 77(^2) = 2, . . . ,it{jn-k) = n — k, 
and also %(ix) = n — k + 1, 7r(i 2 ) = n — k + 2, . . . , 7r(ifc) = n. Therefore the unbounded line-solitons 
of the contour plot of the soliton solution are labeled as [it, n — k + I] for I = 1, . . . , k and [m,j m ] for 
m = 1, . . . , n — k. 

Claim 1. We claim that if . . . , R n are the index sets of the dominant exponentials in the un- 
bounded regions as in Figure [7J then Rg contains I. Therefore Ri+x is obtained from Ri by removing 
I and adding one more index not already in Rg . By Remark 6.4 Claim 1 implies (i?i, . . . , R n ) is the 
Grassmann necklace associated to ir, and therefore implies Theorem |6. 2 1 
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Wc first prove Claim 1 for^ < n — k+1. Clearly 1 £ Ri, since 1 is always the position of an excedance 
of a derangement. Suppose by induction that the claim is true up through £ — 1. Then 

Rt = ((((({it,. . . ,i k } U {jr}) \ {1}) U {j 2 }) \ {2}) • • • U { Jei }) \ {4}. 

Suppose that £ ^ Rg. In steps 1 through I — 1, we have only removed the numbers {1, 2, . . . ,t — 1}, 
and so £ ^ {ii, . . . And we have only added the numbers {ji, . . . , jg~i}, and so £ ^ {jx, . . . , jg-i}. 
Since ^? ^ {ii, . . . we have ir(£) < £, and so £ £ {ji, . . . , Since ^ ^ {ji, . . . , we have 

^ € {jf, j£+i, • ■ • , jn-fe}. But 1 < ji < j 2 < • • • < jn-k and so each element in . . -,j n -k} is 

greater than I. This is a contradiction. 

Claim 2. R n -k+\ = {7r(ii), . . . , ir(ik)}- Note that since Claim 1 is true for i < n — k + 1, R n -k+i 
contains an index for each excedance position i r such that 7r _1 (i r ) < i r < ir(i r ). (These are the elements 
of R\ that remain in each i?2,-R3, ■ • ■ ^Rn-k+i-) Rn-k+i also contains any nonexcedance position j r 
as long as it is not the case that j r — 7r _1 (j s ) for some s. That is, R n -k+i contains any j r such that 
ir~ 1 (j r ) < j r . Therefore we see that R n ^k+i is equal to the set of values that ir takes at the excedance 
positions of 7r. This proves Claim 2. 

We now prove Claim 1 for £ > n—k+1. Again we use induction on £. The claim is true for n — k + 1. 
Suppose that £ ^ Rg but Claim 1 is true for smaller £. Certainly £ £ R n _k+i- So £ ^ Rg means that 
£ must have been removed at some earlier step - say step r, for n — k + 1 < r < £. But the numbers 
removed at these steps were precisely the numbers n—k + l,n — k + 2, ...,£ — 1. This is a contradiction. 
This finishes the proof of Claim 1 and hence of Theorem |6.2| □ 

7. SOLITON GRAPHS ARE GENERALIZED PLABIC GRAPHS 

In this section we will show that we can think of soliton graphs as generalized plabic graphs. More 
precisely, we will associate a generalized plabic graph Pl(C) to each soliton graph C. We then show 
that from Pl{C) - whose only labels are on the boundary vertices - we can recover the labels of the 
line-solitons and dominant exponentials of C. The upshot is that all edge and region labels of a soliton 
graph C may be reconstructed from a labeling of each boundary vertex of C by an integer. 

Definition 7.1. A generalized plabic graph is an undirected graph G drawn inside a disk with n boundary 
vertices labeled 1, . . . , n placed in any order around the boundary of the disk, such that each boundary 
vertex i is incident to a single edge. Each internal vertex is colored black or white, and edges are allowed 
to cross each other in an X-crossing (which is not considered to be a vertex). 

Definition 7.2. Fix an irreducible cell S^ nn of (Gri^n)>o- To each soliton graph C coming from a point 
of that cell we associate a generalized plabic graph Pl(C) by: 

• labeling the boundary vertex incident to the edge {i, iTi} by ~Ki = 7r(i); 

• forgetting the labels of all edges and regions. 

See Figure [8] for a soliton graph C (the same one from Figure [5]) together with the corresponding 
generalized plabic graph Pl(C). 

Remark 7.3. When ir indexes a TP Schubert cell in (Grk. n )>o, the boundary vertices will be labeled by 
1, 2, . . . , n in counterclockwise order, with 1, 2, . . . , n — k labeling the boundary vertices corresponding 
to the y <§; part of the soliton graph. 

We now generalize the notion of trip from [231 Section 13]. 

Definition 7.4. Given a generalized plabic graph G, the trip Ti is the directed path which starts at the 
boundary vertex i, and follows the "rules of the road" : it turns right at a black vertex, left at a white 
vertex, and goes straight through the X-crossings. Note that will also end at a boundary vertex. 
The trip permutation ttq is the permutation such that ^(i) = j whenever Tj ends at j. 
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Figure 8. A soliton graph C and generalized plabic graph Pl(C) for ir = (7, 4, 2, 9, 1, 3, 8, 6, 5). 

We use the trips to label the edges and regions of each generalized plabic graph. 

Definition 7.5. Given a generalized plabic graph G with n boundary vertices, start at each boundary 
vertex i and label every edge along trip Ti with i. Such a trip divides the disk containing G into two 
parts: the part to the left of Ti, and the part to the right. Place an i in every region which is to the left 
of Ti. After repeating this procedure for each boundary vertex, each edge will be labeled by up to two 
numbers (between 1 and n), and each region will be labeled by a collection of numbers. Two regions 
separated by an edge labeled by both i and j will have region labels S and (S \ {i}) U {j}. When an 
edge is labeled by two numbers i < j, we write [i, j] on that edge, or {i,j} or {j,i} if we do not wish 
to specify the order of i and j. 

Theorem 7.6. Consider a soliton graph C — G t (u^) coming from a point A of an irreducible positroid 
cell 5* nn . Then the trip permutation associated to Pl{C) is tt, and by labeling edges and regions of 



Pl{C) according to Definition 7.5 we will recover the original labels in C. 



We invite the reader to verify Theorem |7.6| for the graphs in Figure [8j 



Remark 7.7. By Theorem |7.6| we can identify each soliton graph C with its generalized plabic graph 
Pl(C). From now on, we will often ignore the labels of edges and regions of a soliton graph, and simply 
record the labels on boundary vertices. 

In the proof below, we will sometimes refer to the contour plot from which the soliton graph came; 
it is useful to think about whether edges are directed up or down. 

Proof. We begin by analyzing the edge labels around a trivalent vertex in a soliton graph. They must 
have edge labels [i, j] , [i, m] , and [j, m] in some order, where without loss of generality i < j < m. Recall 
that the slope of a line-soliton labeled [i, j] is n t + Kj. Also recall that we fixed Ki < k 2 < • • • < K n . 
Therefore wc know that the slopes of these three line-solitons are ordered by Ki + Kj < Ki + n m < Kj+K m . 
It follows that a trivalent vertex in the contour plot with a unique edge directed down (respectively, 
up) from the vertex must have line-solitons labeled as in the left (respectively, right) of Figure [9] 

We now fix r between 1 and n, and analyze the set of all edges in the soliton graph C whose label 
contains an r. We aim to show that this set of edges is a trip. 

If r is an excedance value of tt, then we know from Theorem |5 . 5| that there is an edge incident to the 
boundary of C which is labeled [£, r], where £ < r. This is an unbounded edge going to y — > oo in the 
contour plot. And if r is a nonexcedance value, there is an edge incident to the boundary of C which 
is labeled [r, £] where r < £. This is an unbounded edge going to y — > — oo in the contour plot. 

Considering Figure [9] and Definition |7.2| it is clear that the set of all edges containing an r in C will 
be a path between boundary vertices r and 7ry in Pl(C). We call this the soliton path. 
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Figure 9. Contour plots of resonant interactions of three line-solitons 



We now claim that if we start at vertex r and follow the soliton path to vertex n ri then the path 
will have the following property: the path travels down along an edge with labels q and r if and only if 
q < r, and the path travels up along an edge with labels q and r if and only if q > r. 

This claim is clearly true for the first edge of each soliton path. Now we just need to check that the 
claim remains true as we pass through black and white vertices. 

Suppose that we are traveling down along an edge with labels i and r where i < r, and we get to a 
white vertex. Then, looking at the right side of Figure [9j we must have m = r, so the next edge that 
we traverse must be the edge \j, m] in the figure (that is, [j, t]). Note that we will continue to go down 
along an edge with labels j and r, with j < r. 

Suppose that we are traveling down along an edge with labels q and r where q < r, and we get to 
a black vertex. Then, looking at the left side of Figure [9] there are two possibilities. Either we are 
traveling down along the left edge (labeled [j, to] in the figure, so that j — q and m = r), or we are 
traveling down along the right edge (labeled in the figure, so that i — q and j = r). In the first 
case, the next edge we traverse will be the edge labeled [i, m] in the figure, i.e. [i, r], so we will continue 
to go down along an edge with labels i and r, with i < r. In the second case, the next edge we traverse 
will be the edge labeled [j, to] in the figure, i.e. [r, m]. So in this case, our next edge in the path will go 
up along an edge with labels [r, to], where m> r. In all cases, the claim continues to hold. 

There are also three cases to analyze if we go up along an edge. These three cases are completely 
analogous. Therefore the claim is true by induction. 

Finally we note that in all of the above cases, every sequence of edges in the soliton path obeys the 
"rules of the road" . This shows that the soliton paths agree with the trips, completing the proof of 
Theorem 17.61 □ 



8. Contour plots and soliton graphs for positroid cells when t — > ±oo 

Choose a positroid cell Sffi 1 . In this section we will explicitly compute the limit of the contour plot 
C t (uA) (for any A € S'a/T) as * tends to oo or — oo. In particular, we will provide an algorithm that 
constructs the associated soliton graphs, and we will give coordinates for all the trivalent vertices. Most 
of this section will be devoted to the case when t — > — oo, and then we will explain how the same ideas 
can be applied to the case t — > oo. Note that in this section we will assume that x, y, and t are on a 
large scale, while £4, . . . , t m are finite (so we may as well set them equal to 0). 



.1. The definition of the contour plots C± OQ {M.) as t ^ ±00. Recall from Definition 4.1 that 

<tnn 
M 



given A £ Sfy n an d a fixed t £ R, the contour plot C t {uA) is the locus in R 2 where the equation 



k 

,2 1 „3 



f A (x,y) = max <^ ln(Aj(A)Kj) + £ {n h x + K 2 n y + kU) 

JEM i=l 

is not linear. Here we write J = {ji, . . . ,jk}- If t tends to ±00, then each of the terms \n(Aj(A)Kj) 
is negligible. So we just need to analyze 



fM(x,y) = max { 2 (/% 4 z + >c|y + 
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Remark 8.1. The constant terms \a(Aj(A)Kj) provide the information on the locations of line-solitons 



(see (4.2 1), and in particular are responsible for the phase shifts appearing in X-crossings in the contour 
plot Ct{uA) (see [B] for the details). Therefore there will not be any phase shifts in the contour plot 
obtained from Jm{x, y), and the line-soliton of [i, j]-type is simply given by 



x + (m + Kj)y + (nj + Kj + KiKj)t = 0. 



To understand how fM( x jy) behaves as t — > ±oo, let us rescale everything by t. Define x — f and 
y = |, and set 

<t>i{x,y) = Kix + nfy + k - , 

that is, KiX + K?y + aft = t(f>i(x, y). Note that if t is positive then x and y have the same signs as x and 
y, respectively. However, if t is negative then x and y have the opposite signs of x and y. This leads to 
the following definition of the contour plots at t — > oo and t — > — oo. 

Definition 8.2. We define the contour plots C 00 (A4), and C-oo(AA), to be the locus in K 2 where 

max | £ (x, y) | , and mm | ^ 4> h (x , y) | , 
respectively, are not linear. 

Remark 8.3. C 00 (A / i) is the limit of Ct{uA) as t — > oo with the coordinates (x = f,y = f), for any 
A e Sj!% 1 . Similarly, after a 180° rotation, C_ 00 (7W) is the limit of C t (uA) as t -> -oo, for any A e 5^". 
Note that the rotation is required because the positive x-axis (respectively, y-axis) corresponds to the 
negative x-axis (respectively, y-axis). 

Definition 8.4. For 1 < i < j < n, let £y be the line in the iy-plane where </>j(x, y) = <j)j(x, y). And let 
Ui,<,7n be the point where 0,(x, y) = ^(i, y) = m (i, y). 

The following lemma is easy to check. 

Lemma 8.5. Lij has the equation 

x + («j + Kj)y + (k 2 + /tjKj + k 2 ) = 0, 

and the points Vi^ m have coordinates 

Some of the points Vi^ m £ M 2 will be trivalent vertices in the contour plots we construct; such a 
point corresponds to the resonant interaction of three line-solitons of types [i,£], [£,m] and [i,m] (see 
Theorem |8.8| below) . 

8.2. Main results on C± OCl (M.) and their soliton graphs. Consider a positroid cell S^" = S^" 1 
where L is the J-diagram indexing the cell. We will explain how to use L to construct a generalized 
plabic graph G-(L). 

Algorithm 8.6. From a J-diagram L to the graph G_(L): 

(1) Start with a J-diagram L contained in a k x (n — k) rectangle, and use the construction of 
Definition |2. 13| to replace 0's and +'s by crosses and elbows, and to label its border. 

(2) Add an edge, and one white and one black vertex to each elbow, as shown in the upper right 



of Figure [10] Forget the labels of the southeast border. If there is an endpoint of a pipe on 
the east or south border whose pipe starts by going straight, then erase the straight portion 
preceding the first elbow. 
(3) Forget any degree 2 vertices, and forget any edges of the graph which end at the southeast 
border of the diagram. Denote the resulting graph G_(L). 
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(4) After embedding the graph in a disk with n boundary vertices (this is just a cosmetic change 
which we sometimes omit), we obtain a generalized plabic graph, which we also denote G_(L). 
If desired, stretch and rotate G-(L) so that the boundary vertices at the west side of the 
diagram are at the north instead. 

Figure [10] illustrates the steps of Algorithm |8.6| starting from the J-diagram of the positroid cell 
gtnn wnere Tf — (7,4,2,9,1,3,8,6,5). After labeling the edges according to the rules of the road, we 
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Figure 10. Construction of G_(L) where tt(L) 
figure is L. 



(7,4,2,9,1,3,8,6,5). The top left 



will produce the graph from Figure [5] 

Remark 8.7. If every box of L contains a + (that is, 6>£"™ is a TP Schubert cell), then G_(L) will not 
contain any A-crossings. 

The following is the main result of this section. The proof will be given in the next subsection. 

Theorem 



8.6 



to obtain G-(L). Then 



Choose a positroid cell = Sf n = S 1 ™. Use Algorithm 
G-{L) has trip permutation n, and we can use it to explicitly construct C-oo(M) as follows. Label the 
edges of G_(L) according to the rules of the road. Label each trivalent vertex incident to solitons [i,£} 7 
[i,m], and [£,m] by Xi.£ im and give that point the coordinates (x,y) — (KiKt + KiK m + K£K m , — (k, + 



Kit + K m ))- Place each unbounded line-soliton of type so that it has slope «j 



line-soliton of type will automatically have slope Kj 
180° rotation, C_ 00 (A / 1) is the limit of C t (u a) as t - 



.) Recall from Remark 8.3 



oo, for any A g Sffl 1 . 



(Each bounded 
that after a 



Remark 8.9. Although Theorem 8.8 dictates which collections of line-solitons meet at a trivalent vertex, 
it does not determine which pairs of line-solitons form an A-crossing. Which line-solitons form an X- 
crossing is determined by the parameters (ki, . . . , K n ). See Figure 11 for three contour plots based on 
three different choices of (ki, . . . ,K n ). All of them can be constructed using the graph G_(L) from 
Figure [lOj together with Theorem |8.8| 



We can use a very similar algorithm to construct G + (L) from the "dual" J-diagram of L. 



20 



YUJI KODAMA AND LAUREN WILLIAMS 






Figure 11. Three different contour plots C-oc(M) where ir(M) = 
(7,4,2,9,1,3,8,6,5), based on different choices for (k±,...,kq). The left panel 
corresponds to the contour plot with (k±, . . . , Kg) = (—4, —3, —2, —1, 0, 1, 2, 3, 4); note 



that this is the same as G-{L) (cf. Figure 10 1 



Definition 8.10. Given A4 C ('^), we define its dual A4* to be the collection 

M* ={{n+l-j 1 ,n+l-j 2 ,...,n + l-j k } | . . . ,j k } € M}. 

Given it € S n , we define its dual to be the permutation ir* — l o 7r _1 , where i is the involution in S n 
such that l(J) = n + 1 — j. 

Given a J-diagram L, we define its dual to be the J-diagram L* such that 7r(X*) = ir(L)*. 

Remark 8.11. Note that which positroid cell SjK" = S*"" = S£ n ™ a fixed element of (Gr k , n )>o lies 
in depends on a choice of ordered basis (ei,e2, . . . ,e„) for K™. If we relabel each basis element et by 
e-n+i-i, then M., 7r, and L are replaced by their duals M* , ir* , and L*. 
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Figure 12. Construction of G+(L) where tt(L) = (7,4,2,9,1,3,8,6,5). The top left 
figure shows L*. If we compare the unbounded solitons in G+{L) and G-(L), they 
appear to be in a different order. However, their order will be the same at \y\ 3> 0. 
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Theorem 8.12. Choose a positroid cell Sffi = Sf n = Sf n where tt £ S n . Apply Algorithm [sTe] to L* , 
but replace every label j around the boundary of the J-diagram and plabic graph with 7r(n + 1 — j). 
This produces a graph we call G+(L). Then we can explicitly construct Coo(-M) from G+(L), just as 
in Theorem 



See Figure 12 



8.3. The proof of Theorem |8.8[ In this section we present the proof of Theorem |8.8| The main 
strategy is to use induction on the number of rows in the J-diagram L. More specifically, let L' denote 
the J-diagram L with its top row removed. In Lemma 8.14 we will explain that G-(L') can be seen 



.8 correctly 
j(A4), including 



as a labeled subgraph of G-(L). In Theorem 8.17 we will explain that if M! = M.(L'), then there 
is a polyhedral subset of C_ 0o (A / l) which coincides with C- 00 (M I ). And moreover, every vertex of 
C—oo(M) appears as a vertex of C_ 00 (A'l). By induction we can assume that Theorem 
computes C- 00 (M'), which in turn provides us with a description of "most" of C_ 
all line-solitons and vertices whose indices do not include 1. On the other hand, Theorem |5.5| gives a 
complete description of the unbounded solitons of both C_ 00 (A4 / ) and C-oo(M.) in terms of ir(L') and 
7r(L). In particular, C_ 00 (^V{) contains one more unbounded soliton at y 3> than does C_ oc (A^), and 
C-oo(M) contains I more unbounded solitons at y <C where I is the difference in length of the first 
two rows. This information together with the resonance property allows us to complete the description 
of C_ 00 (A / i) and match it up with the combinatorics of G-(L). 



Then Algorithm 8.6 produces a 



Lemma 8.13. Let tt = tt(L) be the derangement associated to L. 
generalized plabic graph G—(L) whose trip permutation is tt. 

Proof. It is clear from the construction that G_(L) is a generalized plabic graph. Note that if we 
follow the rules of the road starting from a boundary vertex of G_(L), we will first follow a "pipe" 
northwest (see the top right picture in Figure |10[ ) , and then travel straight across the row or column 
where that pipe ended. This has the same effect as the bijection of Definition |2.13| □ 

We now present a lemma which explains the relationship between G_(L) and G_(L'), where L' is 
the J-diagram L with the top row removed. 

Lemma 8.14. Let L be a J-diagram with k rows and n — k columns, and let G denote the edge-labeled 
plabic graph associated to L via Algorithm 8.6 Form a new J-diagram L' from L by removing the 



top row of L; suppose that £ is the sum of the number of rows and columns in L' . Let G' denote the 
edge-labeled plabic graph associated to L' , but instead of using the labels {1,2, . . . ,£}, use the labels 
{n — I + l,n — I + 2, . . . , n). Let h denote the label of the top row of L. Then G' is obtained from G 
by removing the trip Th starting at h, together with any edges to the right of the trip which have a 
trivalent vertex on Th- 



We omit the proof of Lemma 8.14| it should be clear after the following example. 



Example 8.15 

Figure |8.6| It illustrates the result of Algorithm 



Figure 13 illustrates Lemma 8.14 with the example of tt = (7,4,2,9,1,3,8,6,5) as in 

applied to the chain of J-diagrams obtained by 



successively adding rows from the bottom of the diagram. We suggest that the reader use the rules of the 



road to fill in all edge labels on these (generalized) plabic graphs. The middle part of Figure 13 gives the 
permutation associated to the corresponding J-diagram. Notice the relationship between the excedances 
in these permutations and the labeled line-solitons on the right side of the figure, e.g. the excedances 
(1, 2, 4, 7) and the soliton index [1, 7], [2, 4], [4, 9], [7, 8] in the top figure. It follows immediately from the 
rules of the road that the sequence of (labeled) plabic graphs on the right side of the figure are nested 
within each other. 

Let . . . ,ik} denote the lexicographically minimal element of Ai. (This corresponds to the col- 
lection of pivots for any A € Sjff 1 .) To simplify the notation, we will assume without loss of generality 
that ii = 1. Now set M! = M(L'). We can also describe M' = {J\ {1} | 1 € J andJ e M}. Our next 
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FIGURE 13. Inductive construction of the generalized plabic graph G-(L) of the case 
7T = (7,4,2,9,1,3,8,6,5). 



goal is to explain in Theorem 8.17 the relationship between C_ 00 (A / i) and C-oo(-M'). However, we first 
prove a useful lemma. 

Lemma 8.16. Consider the point v a ,b,c where 1 (jz. {a, 6, c}. Then at this point, we have that 4>x < <j> a = 
4>b = 4>c- It follows that every region R in C-^Ai) incident to the point v a ,f>, c is labeled by a dominant 
exponential Ej such that 1 e J. 

Proof. Recall that <fii(x,y) = KiX + nfy + «?. A calculation shows that <p a {va.b,c) = 4>b{. v a,b,c) = 

4>c(Va,b,c) = K a K b K c , while ^l(w a ,6,c) = Kl(n a K b + K a K c + K b K c ) + n\(-K a - K b - K c ) + . 

Without loss of generality suppose a < b < c, so then k\ < n a < Kb < k c . It follows that (Kb — 

«i)(ac c — Ki) > 0, which implies that Kik^ + k\k c — k\ < ki,k c . Multiplying both sides by (k — Ki), 
which is positive, we get 

2 2 2 3 

K\K^Kb ~\ K\KqKq K-^K\ K^Kb K-^Kq ~\~ K-^ K^KbK^ K\KbK^. 

Therefore 

Kl(K a Kb + K a K c + KbK c ) K^K a + Kb + K c ) + k\ < K a KbK c , 

which implies that <f>i(v a<biC ) < (f> a (v atb ,c) = 4>b{v a ,b,c) = 4>c{v a ,b, c )- □ 

Theorem 8.17. There is an unbounded polyhedral subset 7?. of C_ 00 (A / t) whose boundary is formed by 
line-solitons, such that every region in 1Z is labeled by a dominant exponential Ej for some J containing 
1. In 1Z 7 C_oo(.M) coincides with C- ocl (M'). Moreover, every region of C_ oc (A^ / ) which is incident to 
a trivalent vertex and labeled by Ej> corresponds to a region of C- oc (M) which is labeled by -Ej'ufi}- 
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FIGURE 14. The contour plot C_ 00 (A^ / ) within the contour plot C-oo(M). Note that 
since t is negative here, —y/t ~> is equivalent to y 3> 0. 



Proof. The proof of the first part of the theorem is straightforward. Note that for any value of y, 
there is an x sufficiently large such that 

(f>i(x,y) « fa(x,y) < • • • < <f> n (x,y). 

This proves the existence of the subset 1Z, where every dominant exponential Ej has the property 
that 1 6 J. Therefore the contour plot within 1Z depends only on the information of M 1 , and hence 
coincides with C-oo(M'). (More specifically, the positions of points and line-solitons are identical, and 
each region label is identical to the one from C_ 00 (A^') except that a 1 is added to the index set.) 

What is less obvious is that essentially every region of C_ 00 (A'1 / ) - that is, every region incident to 
at least one trivalent vertex - also corresponds to a region of C-oo{M). For this we need Lemma [8. 1G 



By definition, all points v a ^, c that appear in C_ oc (A^') have the property that 1 ^ {a, 6, c}. The 
three regions R\, R 2 , R3 incident to w a ,fc,c in C-oo(-M') are labeled by E{J\), E(J 2 ), and E(J 3 ). In 
particular, this means that at region R±, J\ is the subset . . . , jk-i} of M! which minimizes the 
value </)j 1 + ■ ■ ■ + 4>j k _ 1 ■ Without loss of generality we can assume that a £ Ji, b € J2, and c £ J3. By 
there is a neighborhood N of i> a ,&, c where <fii < (j> a . It follows that in NDRi, J\ U{jk = 1} 



5.16 



Lemma 

is the subset of A4 that minimizes the value <pj 1 + ■ ■ ■ + <pj k . Therefore the region Ri of C_oo (A4') which 
is labeled by E Jl corresponds to a region of C_ oc (A^) which is labeled by Ej^rn. Similarly for R 2 and 
R3. This completes the proof of the theorem. 



Figure 14 illustrates how the contour plot C_ 00 (A / / / ) sits inside the contour plot C_oo(A^). Recall 
that Ti represents the trip consisting of all line-solitons labeled [1, j] for any j (cf. Figure [5]). 
Theorem |8 . 1 7| immediately implies the following. 

Corollary 8.18. The set of trivalent vertices in C-^M) is equal to the set of trivalent vertices in 
C-oo(M') together with some vertices of the form i>i & c . These vertices are the vertices along the trip 
T\. In particular, every line soliton in C~ 00 (M) which was not present in C- OCl (M') and is not along 
the trip Tj must be unbounded. And every new bounded line-soliton in C_ co (A / t) that did not come 
from a line-soliton in C_ oc (A / J') is of type for some j. 

We can now complete the proof of Theorem |8.8| This proof will repeatedly use the characterization 
of unbounded line-solitons given by Theorem |5.5| 

Proof. Recall that M = M{L) and M 1 = M(L'), where L' is L with the top row removed. By 



Theorem 8.17 we can construct the contour plot C^ oa (M) inductively from the J-diagram L: we start 
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by drawing the contour plot associated with its bottom row, and then consider what happens when we 
add back one row at a time. On the other hand, by Lemma [8. 14] the construction of Algorithm |8.6| can 
also be viewed as an inductive procedure which involves adding one row at a time to the J-diagram. 



Using Lemma 8.13 and Theorem 5.5 we see that Algorithm 8.6 produces a (generalized) plabic graph 
whose labels on unbounded edges agree with the labels of the unbounded line-solitons for the soliton 
graph of any A <E Sj? n . The same is true for A' G S^™". 

Let us now characterize the new vertices and line-solitons which C_ 00 (A / i) contains, but which 
C-oo(M) did not. In particular, we will show that the set of new vertices is precisely the set of 
v i,b,c (where 1 < b < c), such that either c —> b is a nonexcedance of tt = tt(M), or c — > b is a nonex- 
cedance of tt' = tt(M'), but not both. Moreover, if c — > b is a nonexcedance of tt, then V\^, c is white, 
while if c — > b is a nonexcedance of tt', then v\ t c is black. 



By Corollary 8.18 all new vertices have the form V\ (,, c and lie on the trip T\. Additionally, all new 
line-solitons which begin at some point 1^5 c and which are not on the trip T\ must be unbounded. 
Since the points Wi,6, c are trivalent, each one is incident to either an unbounded line-soliton in C_ 00 (A / J), 
which lies in 1Z C , or is incident to a bounded soliton of type which lies in 1Z. (Possibly both are 
true when i = 1). 

If ^i,b.c is incident to a bounded line-soliton which lies in 1Z, that soliton must have been 
unbounded in C- QO (M'), and hence came from a nonexcedance j — ► i in tt'. (All excedances of tt' are 
also excedances in tt.) In particular, i ^ 1, so we can conclude that v\_b,c — v i,i,j- Conversely, if j — > i is 
a nonexcedance of tt' which is not a nonexcedance of tt, then the corresponding unbounded line-soliton 
from C- oc (M / ) becomes a bounded line-soliton in C- 00 {M) which is incident to v\ i j. This 
characterizes the new points i>i,b, c which are incident to a bounded line-soliton contained in 1Z. 

Each other new point Vxj, c w iU t> e incident to either: 

• one unbounded line-soliton of C- oa (M) which lies in 1Z C (plus two bounded line-solitons 
ofTi), or 

• two unbounded line-solitons of C_ 00 (A'l) which lie in TZ C (plus one bounded line-soliton of Ti). 

Either way, it follows that v\^, c is incident to an unbounded line-soliton where i ^ 1, such that 
j — > i is a nonexcedance of 7r but not a nonexcedance of tt'. Therefore v± t b,c — v i,i,j- 

Conversely, each nonexcedance j — > i of tt (respectively, tt') such that 1 < i < j, and such that j — > i 
is not a nonexcedance of tt' (respectively, tt), gives rise to a point v^j of C_ QO (A^). This is simply 
because these line-solitons must have an endpoint in C_ QO (A^) which did not appear in C^ 00 (M'). 

Also note that if v\^ tC is a new vertex such that c — > b is a nonexcedance of tt, then the line-soliton 
[6, c] must go down (towards y <C 0) from i>i,b )C . However, remembering the resonant condition (see 
Figure |9]), and using the fact that 1 < b < c, we see that [6, c] cannot be the only line-soliton going 
down from v\^.c- Therefore Ui,& )C must have two line-solitons going down from it and one line-soliton 
going up from it, so it is a white vertex. 

Similarly, if Ui,6 jC is a new vertex such that c — > b is a nonexcedance of 7r', then the line-soliton [b, c] 
must go up (towards !/>0) from v\ t b,c- By the resonant condition, we see that [b, c] cannot be the only 
line-soliton going up from wi,b lC - Therefore t>i,6 )C must have two line-solitons going up from it and one 
line-soliton going down from it, so it is a black vertex. 

Using the bijection from Definition [2T3J it is straightforward to verify that the above description also 
characterizes the set of new vertices which Algorithm |8.6| associates to the top row of the J-diagram L. 

Finally, let us discuss the order in which the vertices Ui,6,c occur along the trip T% in the contour 
plot. First note that the trip T\ starts at y -C and along each line-soliton it always heads up (towards 
y 3> 0). This follows from the resonance condition - see Figure [9] and take i = 1. Therefore the order 
in which we encounter the vertices v± & iC along the trip is given by the total order on the y-coordinates 
of the vertices, namely Hi + + K c . 
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We now claim that this total order is identical to the total order on the positive integers 1 + b + c 
that is, it does not depend on the choice of fq's, as long as K\ < ■■■ < K n . If we can show this, then 
we will be done, because this is precisely the order in which the new vertices occur along the trip T\ in 
the graph G-{L). 



v77,i ^ 

A Vjj,k 



Figure 15. 



To prove the claim, it is enough to show that among the set of new vertices c , there are not two 
of the form Vi^t and wij.fc where i < j < k < £. To see this, note that the indices b and c of the new 
vertices i>i,6, c can be easily read off from the algorithm in Definition 2.13 c will come from the bottom 



label of the corresponding column, while b will come from the northwest endpoint of the pipe that i>i,6 )C 
lies on. Therefore, if there are two new vertices and vij.fc, then they must come from a pair of 
crossing pipes, as in Figure [l5j Note that the crossing of the pipes must have come from a in the 
J-diagram. From the figure it is clear that the pipe heading north from the crossing must turn west 
at some point, while the pipe heading west from the crossing must turn north at some point. Both of 
these turning points must have come from a + in the J-diagram, but now we see that the J-diagram 
violates the J-property. This is a contradiction, and completes the proof. □ 



.4. The proof of Theorem 8.12, Theorem 8.12 can be seen as a corollary of Theorem 8.8 
Proof. Recall that «!<•••< n n . We define \i = — K n+ i_i. Then Ai < • • ■ < A n . Set y' — —y. Then 



jem 







— max < 

JEM | 


K k X + Kj.y + Kj\ 


= min < 
JeM 1 




= min < 

JEM | 


k 

E ^n+l-j,X - X n+1 _j.y 


= min 

JEM* 


| E A ; •'• A ;."-v;j 


— min 

JEM* 


jgA^ + A^y' + Afl. 



Therefore Coo(M.) is the locus of M 2 where the last equation above is not linear. Comparing this with 
the definition of C- ao (M), we see that Coo{M.) can be constructed from C_ 00 (A / J*), with each label 
j replaced by n + 1 — j, and with an involution replacing y by —y. The effect of the involution is to 
switch the colors of the black and white vertices in the plabic graph, or equivalently, to replace every 
boundary vertex i of the plabic graph by This completes the proof of the theorem. □ 
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Example 8.19. We invite readers to reconstruct the contour plots in Figure [T] The plots correspond 
to the TP Schubert cell S^ 171 with it = (4,5,1,2,6,3). Take the K-parameters as (ki,...,k 6 ) = 
(— 1, —5,0, 1, 1, |). Calculate the trivalent vertices Ujj.fc = (=F( K i K j + i^j^k + KjKfe),±(«;j + Kj + «;&)) 
obtained from the J-diagram and its dual. There are 8 trivalent vertices for both t ^> and ( < 

one obtains the contour plots for Jm^tV) which 



as shown in Figure 16 Then following Theorem 8 



approximate the plots in Figure [T] 
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Figure 16. The Le-diagrams L, L* , and the plabic graphs for ir = (4,5,1,2,6,3). 
Note that there are two X-crossings corresponding to the O's in L* (cf. Figure [I]). 



9. X-CROSSINGS AND VANISHING PLUCKER COORDINATES 

In this section we show that when the contour plot C t (uA) at \t\ 3> has an X-crossing, then the four 
Pliicker coordinates corresponding to the dominant exponentials incident to that X-crossing satisfy a 
"two-term" Pliicker relation. 

Theorem 9.1. Consider some A £ {Grk in )>o, & n d consider the corresponding contour plot C t (u^) when 
t <C or t 0. Choose l<h<i<j<£<n. In the statements below, S is a (fc — 2)-element subset 
of {1,2, ... ,71} which is disjoint from {h,i,j,£}. 

(1) Suppose there is an X-crossing involving line-solitons [h,£] and 

(a) Then if Kh + ki > Ki + Kj, the dominant exponentials around the X-crossing are as in 



Figure 17 (a). Furthermore, if t < 0, A hes (A) = 0, and if t » 0, A ijS (A) = 0. 



(b) Then if Ki + Kj > Kh + k&, the dominant exponentials around the X-crossing are as in 



Figure 17 (b). Furthermore, if t < 0, A tjS {A) = 0, and if t » 0, A MS (A) = 0. 



(2) Suppose there is an X-crossing involving line-solitons [h,i] and [j,£]. Then the dominant ex- 



ponentials around the X-crossing are as in Figure 17 (c). Furthermore, if t 0, Ajgs{A) = 0, 
and if t > 0, A hiS (A) = 0. 
(3) Suppose there is an X-crossing involving line-solitons [h,j] and [i,£]. Then the dominant ex- 



ponentials around the X-crossing are as in Figure 17 (d). Furthermore, if t <g; 0, Ahjs(A) = 0, 
and if i >0, A l£S (A) = 0. 

In particular, in each of the above cases, we get a "two-term" Pliicker relation. 

• In Case (1), we have A us (A)A hjS (A) = A hlS {A)A jes (A). 

• In Case (2), we have A M s(A)A ljS (A) = A hjS (A)A U s(A). 

• In Case (3), we have A hes (A)A ijS (A) = -A hlS {A)A jes (A). 



Remark 9.2. Notice that Case 3 is impossible for contour plots coming from A E (Grk n )>o, since 
Ah£s(A)Aijs(A) — — Af l is(A)Aj£s(A) implies that one of these four Pliicker coordinates must have a 
sign which is different from the other three. 
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[ij] 



(a) 



[h,l] [/,/] 




(b) 



[h,i] M 




(C) 




(d) 



Figure 17. Different types of A-crossings. 



9.1. Proof of Theorem |9.1[ First note that if we are considering the local neighborhood of an X- 
crossing formed by line-solitons on indices h,i,j,£, then we may as well assume that {h,i,j,£} = 
{1,2,3,4} and S = 0. 

Recall that K\ < ■ ■ ■ < K n . Also recall from Section [8] that 

<j>i(x,y) = KiX + K 2 y + «f, 

and the contour plots C 00 (A / J) and C-oo(M) are defined to be the locus in R 2 where 



JEM 



ma ,x S E <t>ji ( x > v) 



and min < ^ 4>j.(x,y) 



JEM 



are not linear. Finally, for any A £ S^j™, Coo(M.) is the limit of Ct{uA) as t — > oo, and after a 180° 
rotation, C~ 00 {M.) is the limit of Ct{uA) as t — > — oo. 



Recall from Definition 8.4 and Lemma 8.5 the definitions of the lines Ly and the points v a bc — 



{_K a Kb ~\~ K a K c ~\~ KbK Cl — (n a + K b + k c )) e R . The following lemma is obvious. 



Lemma 9.3. Define v v abc = 



— (n a + Kb + k c ). Then vf 2 3 > v 



> v 134 > v 23 



We have the following total order on the slopes of the lines for 1 < i < j < 4: 
(a) If Ki + K4 > k 2 + K3 then 

«3 < k 2 + K3 < ^1 + K4 < k 2 + k 4 < K 3 + K4. 



K3 < K\ + K4 < K 2 + K3 < K 2 + K4 < K3 + K 4 . 



«1 + K 2 < Kl 

(b) If Ki + K4 < K2 + K3 then 

Kl + K 2 < Kl 

Proposition 9.4. If ki + K4 > k 2 + K3 then the configuration of lines Ly for 1 < i < j < 4 is as in the 
left of Figure 18 (up to perturbing the Ki's, which perturbs the slopes of lines while keeping the total 
order as shown above). And if ki + K4 < k 2 + K3 then the configuration of lines is as in the right of 
Figure [18] 

Proof. Let x be the point where L14 and L23 meet. If ki + K4 > k 2 + K3, then L24 must intersect 
both ii4 and L23 above x. This follows from the fact that k 2 + K4 > ki + K4 and > ^234 (fr° m 



Lemma 9.3 1. While if Ki + K4 < k 2 + K3, then L 2 4 must intersect both L 14 and L23 below x. This 

In either case, we can now draw L24, and 



> u 



234- 



follows from the fact that k 2 + K4 > k 2 + K3 and vf 24 
so have locations for the points ^124 and t>234- 

Now consider the placement of L13. If ki + m > k 2 + K3 (respectively ki + K4 < k 2 + K3) then 
Ki + K3 < k 2 + K3 (respectively, ki + K3 < ki + K4). And L13 intersects L14 and L23 in ^134 and W123, 

We now have 



which must satisfy v. 



231 



< V 



131 



< V 



121 



> V 



123' 



So L13 must be as shown in Figure 



locations for all four points ifefe, so we can draw in all six lines Li 



18 



□ 
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Figure 18. The lines L^, based on whether K\ + K4 > n 2 + K3 or n\ + K4 < k 2 + K3. 

Now for each region in the two figures, we will compute the total order on {(j>\, 4>2<, 4>s-, 4>i\- If 
4> a {x,y) < 4>b{x,y) < (f> c (x,y) < <f>d(x,y) then we will write abed as shorthand for this order. Also note 
that if y is finite then for x 3> 0, we have 4>i{x,y) < 4> 2 {x,y) < 4>^(x,y) < 4>4(x,y). This allows us to 
compute the total orders on the 4>iS, as shown in Figure [19} 




Figure 19. The total order on the <fii(x, j/)'s, based on whether k\ + K4 > k 2 + K3 or 

Using Figure [19} we can compute the dominant exponentials. To compute the dominant exponentials 
for t (respectively f > 0) in a given region, we consider the region label abed and choose the 
leftmost two indices (respectively rightmost two indices) such that the corresponding Pliicker coordinate 
is nonzero. 

We now prove Theorem |9.1[ 

Proof. Consider Part (la) of the theorem. Suppose that we see an X-crossing in the contour plot 
involving line-solitons of types [1, 4] and [2, 3]. Let us consider the local neighborhood of this X-crossing, 
looking at the left of Figure [19} Note that in all four regions immediately incident to the X-crossing, 
we have that each of <fii and ^4 is less than each of 02 and ^3. So at t <C 0, if Ai4(A) 7^ 0, then 
this X-crossing would not appear in the contour plot (-E14 would be the dominant exponential in a 
neighborhood of the X-crossing). Therefore we must have Ai4(A) = 0. Now recalling that Cj(ma) for 
t is obtained from C_ 00 (A / J) by a 180 degree rotation, we get dominant exponentials in the four 
regions incident to the X-crossing as predicted by the theorem. 

Similarly, at t 3> 0, if A 2 3(A) ^ 0, then this X-crossing would not appear in the contour plot {E 2 z 
would be the dominant exponential in a neighborhood of the X-crossing.) Therefore we must have 
A 23 (A) = 0. 
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Proving Part (lb) of the theorem is precisely analogous, but we look at the right of Figure 19 
Proving Parts (2) and (3) are very similar, and we leave them to the reader. □ 

10. TP Schubert cells, reduced plabic graphs, and cluster algebras 

The most important plabic graphs are those which are reduced [231 Section 12]. Although it is not 
easy to characterize reduced plabic graphs (they are defined to be plabic graphs whose move-equivalence 
class contains no graph to which one can apply a reduction), they are very important because of their 
application to cluster algebras [57] and parameterizations of cells [33] . 

In this section, after recalling definitions, we will state and prove a new characterization of reduced 
plabic graphs. We then use this characterization to prove that generic soliton graphs for TP Schubert 
cells are in fact reduced plabic graphs. As a consequence, we deduce that the set of dominant expo- 
nentials labeling a generic soliton graph for the TP Grassmannian is a cluster for the cluster algebra 
associated to the Grassmannian. We conjecture that the same is true for TP Schubert cells. 

10.1. Reduced plabic graphs. We will always assume that a plabic graph is leafless, i.e. that it has 
no non-boundary leaves, and that it has no isolated components. In order to define reduced, we first 
define some local transformations of plabic graphs. 

(Ml) SQUARE MOVE. If a plabic graph has a square formed by four trivalent vertices whose colors 
alternate, then we can switch the colors of these four vertices. 



Figure 20. Square move 

(M2) UNICOLORED EDGE CONTRACTION /UNCONTRACTION. If a plabic graph contains an 
edge with two vertices of the same color, then we can contract this edge into a single vertex with the 
same color. We can also uncontract a vertex into an edge with vertices of the same color. 



Figure 21. Unicolored edge contraction 

(M3) MIDDLE VERTEX INSERTION/REMOVAL. If a plabic graph contains a vertex of degree 2, 
then we can remove this vertex and glue the incident edges together; on the other hand, we can always 
insert a vertex (of any color) in the middle of any edge. 



Figure 22. Middle vertex insertion/ removal 

(Rl) PARALLEL EDGE REDUCTION. If a network contains two trivalent vertices of different 
colors connected by a pair of parallel edges, then we can remove these vertices and edges, and glue the 
remaining pair of edges together. 

Definition 10.1. [23] Two plabic graphs are called move- equivalent if they can be obtained from each 
other by moves (M1)-(M3). The move-equivalence class of a given plabic graph G is the set of all plabic 
graphs which are move-equivalent to G. A leafless plabic graph without isolated components is called 
reduced if there is no graph in its move-equivalence class to which we can apply (Rl). 
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Figure 23. Parallel edge reduction 



Theorem 10.2. |23[ Theorem 13.4] Two reduced plabic graphs which each have n boundary vertices are 
move-equivalent if and only if they have the same trip permutation. 

10.2. A new characterization of reduced plabic graphs. 

Definition 10.3. We say that a (generalized) plabic graph has the resonance property, if after labeling 
edges via Definition |7.5[ the set E of edges incident to a given vertex has the following property: 

• there exist numbers i\ < «2 < • • • < i m such that when we read the labels of E, we see the 
labels [«i,i2], [*2,«3], ■ ■ ■ i [im-i>im]> [ii,im] appear in counterclockwise order. 



We call this the resonance property by analogy with the resonance of solitons (see Section 4.2 1. 

Remark 10.4. Note that the graphs in Figure [9] satisfy the resonance property. 

Theorem 10.5. A plabic graph is reduced if and only if it has the resonance property^] 

Remark 10.6. In fact, our proof below also proves that if a generalized plabic graph has the resonance 
property, then it is reduced. 



Proof. By Proposition 



10.8 



below, for every positroid cell S l ^ n there is a reduced plabic graph 
satisfying the resonance property, whose trip permutation equals 7r(L). We will show that the moves 
(Ml), (M2), and (M3) preserve the resonance property. This will show that the entire move-equivalence 
class of Gl satisfies the resonance property. 

i,k] / , \[Lk} 



. U,m) 



(k,m) 



|j,k} /U.ml 



ij,k) 



Ji,k) 



Tj,m) 



Figure 24. The edge-labeling of the square move 



By consideration of the rules of the road, the edges of a plabic graph with a local configuration as 
in the left of Figure [20] must be labeled by the pairs of integers as shown in the left of Figure [24] for 
some i, j, k, m. The edge labeling after a square move is shown in the right of Figure [24) But now note 
that if we compare the top left vertex of the left figure with the bottom right vertex of the right figure, 
their edge labels together with the circular order on them concide. Similarly we can match the other 



three vertices of the left figure with the other three vertices of the right figure in Figure 24 Therefore 
the labels around each vertex at the left of Figure [24] satisfy the resonance property if and only if the 
labels around each vertex at the right of Figure [24] satisfy the resonance property. 

Similarly, the edge labels of a plabic graph with a local configuration as in the left of Figure [2T| must 



be as in the left of Figure 25 for some integers h, i, j, k, m. The right of Figure 25 show the new edge 
labels we'd get after a unicolored edge contraction. Note that in order for the configuration on the left 
to satisfy the resonance property, we must have h, i,j, k, m be cyclically ordered, e.g. h<i<j<k<m 



Recall from Definition 



that our convention is to label boundary vertices of a plabic graph 1,2, ... ,n in counter- 



clockwise order. If one chooses the opposite convention, then one must replace the word counterclockwise in Definition 



10.3 by clockwise. 
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Figure 25. The edge-labeling of the unicolored edge contraction 



or i < j < k < m < h or j < k < m < h < i or .... Similarly, in order for the configuration at the right 
of Figure [25] to satisfy the resonance property, we must have h, i, j, k, m be cyclically ordered. Therefore 
the move (M2) preserves the resonance property. 



The move (M3) trivially preserves the resonance property. All edges in Figure 22 will be labeled 
for some i and j. Therefore we have shown that moves (Ml), (M2), and (M3) preserve the resonance 
property. 

By |23l Theorem 13.4], for any two reduced plabic graphs G and G' with the same number of 
boundary vertices, the following claims are equivalent: 

(1) G can be obtained from G' by moves (M1)-(M3) 

(2) G and G' have the same (decorated) trip permutation. 

Therefore it follows that all reduced plabic graphs with the (decorated) trip permutation 7r satisfy the 
resonance property. Letting L vary over all J-diagrams, we see that all reduced plabic graphs satisfy 
the resonance property. 

Now we need to show that if a plabic graph G has the resonance property, then it must be reduced. 
Assume for the sake of contradiction that it is not reduced. Then by [231 Lemma 12.6], there is another 
graph G' in its move-equivalence class to which one can apply a parallel edge reduction (Rl). Since 
applications of (Ml), (M2), (M3) preserve the resonance property, G' has the resonance property. But 
then, it is impossible to apply (Rl). This is because an edge-labeling of a local configuration to which 

[ij] 

Figure 26. The edge-labeling of the parallel edge reduction 



one can apply (Rl) must be as in the left hand side of Figure 26 However, this local configuration 
violates the resonance property: it has a trivalent vertex with two incident edges which have the same 
edge-label. Therefore one cannot apply (Rl) to G' so G must be reduced. □ 

We now provide an algorithm from [231 Section 20] for associating a reduced plabic graph G^ to any 
J-diagram L. The plabic graph will have the trip permutation tt(L). 

Algorithm 10.7. [H Section 20] 

(1) Start with a J-diagram L contained in a k x (n — k) rectangle, and label its southeast border 
from 1 to n, starting from the northeast corner of the rectangle. Reflect the figure over the 
horizontal axis. 

(2) From the center of each box containing a +, drop a "hook" up and to the right, so that the 
arm and leg of the hook extend beyond the east and north boundary of the Young diagram. 
Consider the hook graph H(L) formed by the set of all such hooks. 



(3) Make local modifications to H(L) as in Figure 27 
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Figure 27. 



(4) The border labels from the first picture become labeled "boundary vertices;" they are labeled 
1 to n in counterclockwise order. After embedding the figure in a disk, we have a plabic graph 
which we denote by G^. 

Figure [28] illustrates the steps of Algorithm |10.7| 
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Figure 28. A plabic graph for <S* n ™ with tt = (6, 4, 2, 7, 1, 3, 8, 5) 



Proposition 10.8. The plabic graph G^has the resonance property. 

Proof. We prove this directly by analyzing the trips in the graph. One can collapse the plabic graph 
G^ back to the hook graph H(L), and consider how the trips look in H(L). In H(L), the trips have 
the following form: 

• a trip which starts from the label of a vertical edge in the original J-diagram first goes west as 
far as possible, and then takes a zigzag path north and east, turning whenever possible. 

• a trip which starts from the label of a horizontal edge in the original J-diagram first goes south 
as far as possible, and then takes a zigzag path east and north, turning whenever possible. 



See Figure 29 for a depiction of the general form of the trips, as well as the trip which begins at 8 in 
the example from Figure [28] 

We now analyze the edge labeling around a black vertex in G^ which comes from a trivalent vertex 



in the hook graph H(L), see Figure 30 By consideration of the zigzag shape of the trips, the trip which 
approaches the black vertex from above must come straight south from some boundary vertex labeled 
j, while the trip which approaches the black vertex from the right must come straight west from some 
boundary vertex i. The trip which approaches the black vertex from below could have started from a 
boundary vertex k which is either southeast of i or west of j, see the first two pictures in Figure [30] 
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Figure 29. 



Either way, the resulting edge labeling will be as shown in the third picture in Figure [30} Clearly i < j, 
and either k < i or k > j. Therefore the edge labeling in the third picture in Figure |30| satisfies the 
resonance property. 



-d k 



,r wl 



-< k 



J 






(J.k) 


(i,k) 


, (ij) 



Figure 30. 



The argument for a white vertex in which comes from a trivalent vertex in the hook graph H(L) 
is analogous; see the fourth, fifth, and sixth pictures in Figure [30} 

] _ T 

Finally we analyze the edge labeling around a pair of vertices in which came from a degree 4 
vertex in H(L). In H(L), the trip which approaches the vertex from above must come straight south 
from some boundary vertex labeled j, while the trip which approaches the vertex from the right must 
come straight west from some boundary vertex i. As before, i < j. Let k and £ denote the boundary 
vertices whose trips approach the degree 4 vertex from the left and below, respectively. The resulting 



edge-labeling is shown at the right of Figure 31 



There arc multiple possibilities for the trips starting from k and I; Figure 31 shows several of them. 
The only restriction is that neither k nor I lies in between i and j. In particular, one of the following 
must be true: 



• k < i < j and £ < i < j 

• i < j < k and i < j < £ 

• k < i < j < £ 

• £ < i < j < k 



In all cases, the edge labeling in Figure 31 satisfies the resonance property. □ 



0,1) 



Figure 31. 
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Corollary 10.9. Suppose that A lies in a TP Schubert cell, and that for some time t, Gt(uA) is generic 
with no X-crossings. Then G t {uA) is a reduced plabic graph. 



Proof. By Remark |7.3| and Theorem |7.6( every such graph is a plabic graph which satisfies the 
resonance property. □ 

10.3. The connection to cluster algebras. Cluster algebras are a class of commutative rings, intro- 
duced by Fomin and Zelevinsky , which have a remarkable combinatorial structure. Many coordinate 
rings of homogeneous spaces have a cluster algebra structure: as shown by Scott |27j . the Grassmannian 
is one such example. 

Theorem 10.10. [27] The coordinate ring of (the affine cone over) Grk, n has a cluster algebra structure. 
Moreover, the set of Plucker coordinates whose indices come from the labels of the regions of a reduced 
plabic graph for (Grk, n )>o comprises a cluster for this cluster algebra. 

Remark 10.11. In fact [27J used the combinatorics of alternating strand diagrams, not reduced plabic 
graphs, to describe clusters. However, alternating strand diagrams are easily seen to be in bijection 
with reduced plabic graphs [23] . 

Theorem 10.12. The set of Plucker coordinates labeling regions of a generic trivalent soliton graph for 
the TP Grassmannian is a cluster for the cluster algebra associated to the Grassmannian. 



Proof. This follows from Corollary 10.9 and Theorem |1 0.10 □ 



Conjecturally, every positroid cell <S>£"" of the totally non-negative Grassmannian also carries a cluster 
algebra structure, and the Plucker coordinates labeling the regions of any reduced plabic graph for 5*"™ 
should be a cluster for that cluster algebra. In particular, the TP Schubert cells should carry cluster 
algebra structures. Therefore we conjecture that Theorem |10.12| holds with "Schubert cell" replacing 



"Grassmannian." Finally, there should be a suitable generalization of Theorem 10.12 for arbitrary 
positroid cells. 

11. The inverse problem for soliton graphs 

The inverse problem for soliton solutions of the KP equation is the following: given a time t together 
with the contour plot Ct{uA) of a soliton solution, can one reconstruct the point A of Grk, n which gave 
rise to the solution? Note that solving for A is desirable, because this information would allow us to 
compute the entire past and the entire future of the soliton solution. 

We remark that in this section we work with the contour plots C t (uA), as opposed to their limits 
C ±00 (.M), because the former include the information of the values of the Plucker coordinates. However, 
our description of the limits C±ao{Ai) will be a crucial tool for solving the inverse problem when \t\ is 
sufficiently large. 



In this section we will solve the inverse problem whenever A S (Grk, n )>o (see Theorem 11.2 1, or 



whenever A <E (Grk, n )>o an d \t\ is sufficiently large (see Theorem 11.3) 



Lemma 11.1. Fix ki < ■ ■ ■ < K n as usual. Consider a contour plot Ct(uA) of a soliton solution coming 
from a point A of (Gr£, „)> for any time t. Then from the k^'s, the contour plot, and the time t 
we can identify what cell S l ™ n the element A comes from, and reconstruct the labels of the dominant 
exponentials and the values of all Plucker coordinates corresponding to these dominant exponentials in 
the contour plot. 

Proof. Since the pairwise sums of the k^s are all distinct, we can determine from the contour plot 
precisely how to label each line-soliton with a pair Note here that some of the edges in C((wa) 

correspond to the phase shifts. However, those can be easily identified by checking the types of solitons 
at the intersection point, since a phase shift appears as the interaction of two solitons of [i, j}- and [£, m]- 
types with either i < j < £ < m or i < £ < m < j. From the labels of the unbounded line-solitons, we 
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can use Theorem |5.l| t o determine which positroid cell S% in the element A comes from. We then use 
Definitions |4.8| and |7.2| to pass from the contour plot to the generalized plabic graph. Finally we use the 



rules of the road (Definition 7.5 and Theorem 7.6) to label the dominant exponentials in each region of 
the contour plot. Alternatively, we can label the dominant exponentials by using Lemma |4.5| together 
with the fact that for i<0, the dominant exponential is Ej, where Aj is the lexicographically minimal 
Pliicker coordinate which is nonzero on 5£ nn (/ is the set of excedance positions of 7r). 



Now recall that the equation of each line-soliton is given by (4.2 1. Since each line-soliton has been 
labeled by and we know the re^'s, we can solve for all ratios of Pliicker coordinates which appear 
as labels of adjacent regions of the contour plot. The Pliicker coordinates are only defined up to a 
simultaneous scalar, so without loss of generality we set the Pliicker coordinate Aj with lexicographically 
minimal index set I equal to 1. □ 

Using the cluster algebra structure for Grassmannians, we can now prove the following. 

Theorem 11.2. Consider a generic contour plot C^ua) of a soliton solution which has no A-crossings, 
and which comes from a point A of the totally positive Grassmannian at an arbitrary time t. Then 
from the contour plot together with t we can uniquely reconstruct the point A. 



Proof. By Lemma 11.1 we can compute the labels of the dominant exponentials in the contour plot 



and the values of the corresponding Pliicker coordinates. And by Theorem 10.12 the set of dominant 
exponentials labeling Ct(uA) forms a cluster c for the cluster algebra A associated to the Grassmannian. 
Since the coordinate ring of the Grassmannian is a cluster algebra (whose cluster variables include the 
set of all Pliicker coordinates) , it follows that we can express each Pliicker coordinate of A as a Laurent 
polynomial in the elements of c. Therefore we can determine the element A £ {Grk : n)>o itself. □ 

Theorem 11.3. Fix k\ < ■ ■ ■ < n n as usual. Consider a generic contour plot Ct(iiA) of a soliton solution 
coming from a point A of a positroid cell S^™", for i <§; 0. Then from the contour plot together with t 
we can uniquely reconstruct the point A. 

Theorem |8.8| w ill be instrumental in proving Theorem |11.3| However, we first remind the reader 
(see Remark |8.9[ ) that the combinatorics of the A-crossings in the contour plot may differ from the 
combinatorics of the A-crossings in the graph G_(L). We will describe these differences using the 
following notion of slide. 

Definition 11.4. Consider a generalized plabic graph G with at least one A-crossing. Let ii a ,6,c be a 
trivalent vertex (with edges labeled [a, b], [a, c], and [b, c]) which has a small neighborhood N containing 
one or two A-crossings with a line labeled but no other trivalent vertices or A-crossings. Here 

{a,b,c} and must be disjoint. Then a slide is a local deformation of the graph G which moves 

the line so that it intersects a different set of edges of v a ,b,c, creating or destroying at most one region 
in the process. 

Lemma 11.5. Consider the contour plot Ct{uA) of a soliton solution coming from A £ S^™" at t <C 0. 
Then its soliton graph either coincides with G-(L) or differs from it via a series of slides. Similarly, 
if one considers two contour plots of a soliton solution coming from A £ S 1 ™ at t <C 0, which are 
computed using different sets of parameters . . . , K n ), then one can be obtained from the other via 
a sequence of slides. 



Proof. Theorem 8.8 gives a precise description of the trivalent vertices which appear in the contour 
plot, together with the topology of how they are connected to each other (including the circular order 
of the three line-solitons [a, 6], [a,c], and [6, c] incident to a given vertex u a ,t> )C ). In particular, the 
combinatorics of the trivalent vertices and their incident edges in the contour plot is precisely that 
which appears in G_ (L) . 
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FIGURE 32. Examples of slides. These contour plots correspond to the same Le- 
diagram L with ir(L) = (5, 3, 2, 1, 4), but these plots all differ from G-(L). 



The only feature of the contour plot which Theorem |8.8| does not determine is which pairs of line- 
solitons form an X-crossing. Therefore the topology of the contour plot may differ from that of G_ (L) 
via a sequence of slides. In each slide, a line-soliton of type [i, j] may pass across a trivalent vertex u 0) & )C , 
changing the location of the X-crossing, or possibly replacing one X-crossing with two X-crossings (or 
vice-versa). Note that the indices {i,j} and {a, b, c} must be disjoint, since an X-crossing involves two 
lines-solitons with disjoint indices. □ 

Lemma 11.6. Consider two contour plots which differ by a single slide. Let Si and S2 denote the 
two sets of Pliicker coordinates corresponding to the dominant exponentials in the two contour plots. 
Then from the values of the Pliicker coordinates in Si, one can reconstruct the values of the Pliicker 
coordinates in S2, and vice- versa. 



Proof. Recall from Theorem |9.1| that the four Pliicker coordinates incident to an X-crossing are 
dependent, in particular they satisfy a "two-term" Pliicker relation. Now it is easy to verify the lemma 
by inspection, since each slide only creates or removes one region, and there is a dependence among the 
Pliicker coordinates labeling the dominant exponentials. The reader may wish to check this by looking 
at the first and second, or the second and third, or the third and fourth contour plots in Figure |32} □ 

Proposition 11.7. Consider a generic contour plot C*(ua) of a soliton solution coming from a point A 
of a positroid cell Sffi = Sf n = S£ nn C (Gr k . n )> for t < 0. Define 

T={I\ JG (^j and£; is a dominant exponential in Ct{uA)}- 

If we are given 7r and the values of each Aj(A) for I ET~, then we can determine A. 

Proof. For t -C sufficiently small, the contour plot C t {uA) can be approximated by C_ 00 (M). 
Therefore Ct(uA) will either have the same topological structure as the graph G-(L), or by Lemma 



11.5 it will differ from G—(L) via a sequence of slides. And so by Lemma 11.6 we can determine the 
values of the Pliicker coordinates labeling the regions in G_(L). 

We now use Theorem |11.15| below, which shows that the set of Pliicker coordinates labeling the 



regions in G_ (L) contains the set of Pliicker coordinates described in Theorem 11. 9| (the latter result 



is due to Talaska [29 ). By Theorem 11.9 from this set of Pliicker coordinates, one can reconstruct all 
nonzero Pliicker coordinates. This determines the point A € (Grfe )W )>o- □ 

We now turn to the proof of Theorem 11. 3[ which follows easily from our previous results. 



Proof. By Lemma [ll.l[ we can identify the cell S^™" = S^ 171 containing A. We can also compute 
the labels of the dominant exponentials in the contour plot and the values of the corresponding Pliicker 



coordinates. But now by Proposition 11.7 we can determine the element A g (Gr^ „)> . □ 
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Since we also have the generic contour plot for t ^> 0, we can solve the inverse problem for this case 
as well, using now the dual J-diagram. 

Theorem 11.8. From a generic contour plot C t {uA) for t 3> 0, we can uniquely reconstruct the point A. 

Proof. This is equivalent to Theorem |11.3[ since the dual J-diagram is a relabeled J-diagram. □ 

11.1. Reconstructing a point of (Gr^ ,„)>o from a minimal set of Pliicker coordinates. Talaska 
[2"9"] studied the problem of how to reconstruct an element A 6 (Grfc i?l )> from a subset of its Pliicker 
coordinates A/ (A). For each cell Sjtfp, she characterized a minimal set of Pliicker coordinates which 
suffice to reconstruct the corresponding element of Sjpj™. Her proof worked by explicitly inverting 
Postnikov's boundary measurement map [23]. In this section we review some of Talaska's work. 

Given a J-diagram L of shape A which fits inside a k x (n — k) rectangle, we construct a planar 
network Nl as follows. First we draw a disk whose boundary consists of the north and west edges of 
the k x (n — k) box and the path determining the southeast boundary of A. Place a vertex, called a 
boundary source, at the end of each row, and a vertex, called a boundary sink, at the end of each column 
of A. Label these in sequence with the integers {1,2, ... ,n}, following the path from the northeast 
corner to the southwest corner which determines A. Let / = {ii < i 2 < • • • < ik} be the set of boundary 
sources, so that [n] \ I = {ji < j% < ■ ■ ■ < j n -k} is the set of boundary sinks. 

Given a box b in L which contains a +, we drop a hook down and to the right, extending the two 
segments comprising the hook all the way to the boundary of the disk. We direct the horizontal segment 
left and the vertical segment down. After forgetting the +'s and 0's in L, we now have a planar directed 
network. There is exactly one face for each box 6 in A containing a +, and in addition, there is one face 
whose northwest boundary is the boundary of the disk. See the first two pictures in Figure |33| which 
shows this construction applied to the J-diagram from Figure [TO] 
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Figure 33. A minimal set of Pliicker coordinates for n = (7, 4, 2, 9, 1, 3, 8, 6, 5). 

Talaska proved the following. 

Theorem 11.9. [HI Corollary 4.2], [5H1 Lemma 2.1] Let L be the J-diagram of a positroid cell S^™", 
and A e 5^ nn . Define a set of Pliicker coordinates T(L) := {Aj(A)} U {Aj^(A)}b, where b ranges 
over all boxes containing a + in the J-diagram, and J(b) is the destination set of the northwest-most 
path collection lying weakly southeast of the b hook. Then T(L) is a totally positive base; that is, each 
nonzero Pliicker coordinate on S 1 ™ 1 can be written as a subtraction-free expression in the elements of 
T(L). It follows that one can reconstruct A from T(L). 



The right side of Figure 33 shows the Pliicker coordinate Ajym associated with each box b which 



contained a + in the J-diagram. 

Definition 11.10. Given a J-diagram L and a box b containing a +, let if, be the J-diagram obtained 
from L by changing the content of each box to a if the box lies in a row above & or a column left of b. 
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Definition 11.11. Wc define a total order on elements of (^) by saying that J 1 = {j\ > • • • > j\} > 
J 2 = {j 2 > • • • > j 2 } if in the first position i where they differ, jj > j 2 . Then if M C (^), we refer to 
the largest element of M. as lexicographically maximal. 

The following lemma is implicit in |29j . 



Lemma 11.12. Let Ai :— M.(Lb) be the collection of Pliicker coordinates which are positive on the cell 
Then Aj(m is the lexicographically maximal element of M.. 

The following lemma is obvious. 



Qtnn 
Lb 



Lemma 11.13. The generalized plabic graph G-(Lb) that Algorithm 
inside G-(L). 



6 associates to Lh is contained 



Lemma 11.14. Consider the graph G-(L) constructed by Algorithm 8.6 Let J be the fc-element subset 
of [n] which labels the region R that comes from the northwest corner of the J-diagram L (see the 
bottom left picture in Figure 10). Then J is the lexicographically maximal element of A4(L). 



Proof. One may prove directly that the label of R coincides with the destination set of the northwest- 
most path collection in Nl. 

Alternatively, we may prove this using a soliton argument. Let A E S L nn . When x 0, we have 
that E n > E n -i > ■ ■ ■ > E\. Therefore if J is the lexicographically maximal element of M.(L), then 
the term Aj(A)Ej dominates the T-function ta- It follows that Ej is the dominant exponential of any 
contour plot Ct(uA) in the region where x 3> 0. By Lemma 11.5 for t <C 0, this contour plot coincides 
with G_(L) up to a series of slides, none of which will change the label of the region at x >> 0. And 
this region corresponds to the region coming from the northwest corner of the J-diagram L. □ 

Theorem 11.15. Let L be a J-diagram. Consider the set S of Pliicker coordinates associated to the 



dominant exponentials in the regions of the graph G_(L), as constructed in Algorithm 8.6 Then S 
contains T(L). In particular, one may reconstruct the element A £ S L nn from the values of the Pliicker 
coordinates labeling G_(L) for t <C 0. 



Proof. By Lemmas 



11.12 



11.13 and 



11.14 



for each box b in L, Ej^) is a dominant exponential 
labeling a region of G-(L). Also, Wl is the lexicographically minimal element of Ai(L), then Ej labels 
the region at x <C in G_(L). Therefore T(L) C S. The proof now follows from 11.9 □ 



12. Soliton graphs for (Gr 2i „)>o and triangulations of an b-gon 

The main result of this section is an algorithm for producing all soliton graphs for (Gr 2i „)>o- We 
state the algorithm and main theorem in Section |12.1[ and then prove it in the rest of the section. 



12.1. Triangulations of an n-gon and soliton graphs for (Gr2, n ) 



>o- 



Algorithm 12.1. Let T be a triangulation of an n-gon P, whose n vertices are labeled by the numbers 
1,2, ... ,n, in counterclockwise order. Therefore each edge of P and each diagonal of T is specified by 
a pair of distinct integers between 1 and n. The following procedure yields a labeled graph ^(T). 

(1) Put a black vertex in the interior of each triangle in T. 

(2) Put a white vertex at each of the n vertices of P which is incident to a diagonal of T; put a 
black vertex at the remaining vertices of P. 

(3) Connect each vertex which is inside a triangle of T to the three vertices of that triangle. 

(4) Erase the edges of T, and contract every pair of adjacent vertices which have the same color. 
This produces a new graph G with n boundary vertices, in bijection with the vertices of the 
original n-gon P. 
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(5) Add one unbounded ray to each of the boundary vertices of G, so as to produce a new (planar) 
graph ty(T). Note that \&(T) divides the plane into regions; the bounded regions correspond to 
the diagonals of T, and the unbounded regions correspond to the edges of P. 



See Figure 34 Our main result in this section is the following. 



Theorem 12.1. The graphs 9(T) constructed above are soliton graphs for {Gr2 >n )>0i an d conversely, 



up to (M2)- equivalence (see Section 10.1) any trivalent generic soliton graph for (Gr2, n )>o comes from 



this construction. Moreover, one can realize each graph vf'(T') by either: 

• choosing an arbitrary A 6 (Gr2. n )>a o,nd varying the higher times t^ 



• fixing an arbitrary collection of higher times ts, 
appropriate A G (G , r 2 . Il )>o- 



, t n appropriately, or 
,t n , and using the torus action to choose an 





Soliton graph 



Contour plot 



FIGURE 34. Algorithm 12.1[ starting from a triangulation of a hexagon. The right 
figure shows the corresponding contour plot for the soliton solution for (GV 2j 6)>o- 



Remark 12.2. The process of flipping a diagonal in the triangulation corresponds to a mutation in the 
cluster algebra. In the terminology of reduced plabic graphs, a mutation corresponds to the square 



move (Ml) (see Section 10.1). In the setting of KP solitons, each mutation may be considered as an 



evolution along a particular flow of the KP hierarchy defined by the symmetries of the KP equation. 

Remark 12.3. It is known already that the set of reduced plabic graphs for the TP part of Gr 2 . n all 
have the form given by Algorithm 1.1. And we have proved a theorem which says that every generic 
soliton graph is a reduced plabic graph. Therefore from that theorem it follows immediately that every 
soliton graph for the TP part of Gr2, n must have the form of Algorithm 1.1. However, in the rest of 
this section, we will not rely on that theorem, and will instead prove Theorem |12.1| directly. 

12.2. Soliton graphs and flags of exponentials. Recall that the r-function for the TP cell of Gr2, n 
is given by 

r A = £ A iJ (A)E i>j , 

l<i<j<n 

where Eij = («j — K^E^Ej, Ej = exp9j, and 

n 

0j = ^2 Kj t m . 

m—1 
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Definition 12.4. Given a contour plot Ct {uA) of a soliton solution, we associate a flag of exponentials 
I3 C I4 C • • • C I n = [n] where Ij = . . . , ij} c [n] as follows. Choose a trivalent vertex of C tQ {uA), 
and let £'i 1 ,i 2 , Ei 2i i 3 , and Ei 1 ^ 3 denote the three dominant exponentials labeling the regions around the 
vertex. Set ^3 = {11,12,23}. Then find a region adjacent to one of these three. Its label has the form 
EijM (f° r 1 < i < 3), so set I4 = {ii, i%, 13, 14}. At the fcth step, find a region R adjacent to the set of 
regions whose dominant exponentials involve only {Ei t , . . . , Ei k l }. The label of R has the form Ei jtm 
where 1 < j < k — 1, and m ^ {ii, j 2 , . . . , ik-i}- Let i^ — m, and set = {i±, i 2 , . . . , ik}- 

Remark 12.5. Note that one contour plot will in general give rise to many different flags of exponentials. 
However, we will show that one can reconstruct the soliton graph of C t[) (uA) from any one of its flags 
of exponentials. 

Lemma 12.2. Let X = J3 C I4 C • • • C /„ = [n] be any sequence of subsets of [n] where |Jfc| = k. Write 
Ik = {iii *2> • • • ) *fc}- We have the following: 

(1) Fix A £ (G' , 2,ra)>o- Then by varying the higher times t 3 , £4, . . . ,t n we can realize I as a flag of 
exponentials for a contour plot for A. 

(2) Fix the higher times £3, £4, . . . , i„. Then by using the torus action on (Gr2, n )>o, we can realize 
I as a flag of exponentials for a contour plot for some A £ (Gr2, n )>o- 

Proof. Note that by varying the higher times tj, j = 3, 4, . . . , n, we can adjust each 6j (used in the 

n 

definition of Ej = exp#j with 9j = ^ K™i m ) so that for some finite region in the icy-plane, we have 

m— 1 

9 H w 6 i2 « 9 l3 > 6>, 4 > ■■• > in . 

It follows that 

(12.1) E n a £ l2 a £ i3 > >•••> £ in . 

On the other hand, if we keep the higher times fixed, but scale the columns of A £ (Gr2, n )>o by a torus 
element (li,£2, ■ ■ ■ ,£ n ) such that 

~ ~ ^13 ^24 ' ' ' 5 



then we may also assume that (12.1 ) is satisfied. Namely, we can scale the elements Ej instead of the 
columns of A, writing 8j = KjX + Kjy + 9® where tj = exp9j. 

It follows from ( 12.1 ) that for each I > 3 we can approximate the r- function by 

ta « A := £ &*j(A)EiJ, 
i<jeii 

and increasing £ gives a better and better approximation. 

Note that the first approximation t\ corresponds to a local region in Ct {uA) around a trivalent 
vertex, and each subsequent approximation corresponds to increasing that region by adding the contri- 
bution from one more exponential. Therefore X is a flag of exponentials for C to (u^). □ 

From a flag of exponentials for a contour plot, we will show that we can construct its soliton graph. 
Note that when we approximate ta by t a — £ Aij(A)Eij, the corresponding contour plot for t a is 

the contour plot for the r-function coming from the matrix obtained from A by forgetting all columns 
whose indices are outside the set Ig. In particular, this is a contour plot for {Gr 2 /)>o, which appears 
inside the original contour plot Ct (uA)- 

To reconstruct the soliton graph from a flag X, we start with the following lemma. 

Lemma 12.1. Every soliton graph for {Gr2 t n)>o contains a trivalent vertex whose edges are line-solitons 
of types [1, j], [j,n], [l,n] for some 1 < j < n. 
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Proof. Recall that by Theorem 5.5 every soliton graph for (Gr2, n )>o has two unbounded line- 
solitons of types [1, n — 1] and [2, n] (in the j/> direction), and n — 2 unbounded line-solitons of types 
[1,3], [2,4],..., [n- 2, n]. 

Consider a contour plot Ct (uA) for (Gr2 t n)>o- We can choose for it a flag of exponentials I = I3 C 
■ ■ • C in such that there is a k such that 1 € Ik but n £ Ik- 

Choose I such that 1 € It and It+i = Ii U {n}. Order the elements of Z i so that we have 7^ = {1 < 
•••<&< c}. Then 



5.1 



the contour plot for r 



^4 has 



{1 < • • • < 6 < c < n}. Then by Theorem 
two unbounded line-solitons at y S> 0, of which the rightmost one is an unbounded line-soliton of type 
[1,6]. That line-soliton has the region Eb c to its right. The contour plot also has £ — 2 line-solitons at 
i/CO. The contour plot for t^ 1 has the same I — 2 unbounded line-solitons at y <C 0, plus one more 
unbounded line-soliton of type [b,n] at the right. (The [b,n] line-soliton has the region E^ c to its left.) 
On the other hand, the contour plot's two unbounded line-solitons at y 3> are different. In particular, 
the one on the right is of type [1, c]. Note that the bounded portion of the contour plot for t a must 
appear inside the contour plot for t a +1 . 

Let us analyze how the contour plot grows when we pass from t a to t a +1 . Our goal is to prove: 

Claim 12.3. In the contour plot for t a +1 at y ^> 0, a soliton of type [1,6] intersects a soliton of type 



If we know the claim is true, then we will be done, because an intersection of solitons of types [1, 6] 
and [6, n] will produce a third soliton of type [1, n] by the resonance property. 

To prove the claim, let v be the vertex at the top of [1, 6] in the soliton graph. Note that the [1, 6] 
soliton is next to the region E^ c , Without loss of generality v is trivalent, so at v, the [1, 6] soliton meets 
two solitons of types {l,x} and {a;, 6}. 




[Ml 




Figure 35. Cases 1 and 2. 
There are two cases, either x > b or x < b. In Case 1 (resp. Case 2), by the resonance property, the 



local configuration in the soliton graph must look like the left side (resp. right side) of Figure 35 



In the first case, x cannot be c, since [x, 6] is next to region Eb^ c . Therefore x = n and we are done. 
In the second case, [l,x] must be next to region Ej, tC , so one index in {l,x} must coincide with one 
index in {6, c}. But this is a contradiction, since x < b means that x cannot be equal to either 6 or c. □ 

Corollary 12.1. Given a contour plot for {Gr2^ n )>o, we can always choose a flag of exponentials X 
such that the initial index set is I3 — {1, j,n} for some j. 

Proposition 12.1. LetX be a flag of exponentials for a contour plot Ct (u a)- Order the indices in the 
set Ii = {1, j 2 , j 3 , . . . , ji-\,n) by 



Write Ii + i 



1 = h < 32 < 33 < 
Ii U {m}, where m satisfies 

1 = ii < h < ■ ■ • < is < m < j s+ i < 



< ji-i < 3i 



< Jl-i < Jl 
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Then when we go from the soliton graph for to the soliton graph for t_^ +1 7 we add precisely one 
bounded region whose dominant exponential is Ej 3 j . That is, adding one column with index m to 
the A-matrix has the effect of adding to the previous graph one bounded region with label Ej s j , where 
js < m < j s +i- We also add two new bounded line-solitons of types [m,j s +i] and [j Sl m] incident to 
Ej^j g+1 , plus three unbounded line-solitons of types [j' s _i,m], [j s ,j s +i], [to, js+2], o,s in the configuration 
shown in Figure \30^ 

Proof. Since the submatrix of A with column set In is in (Gr 2 .^)>o, the unbounded solitons in the 
soliton graph for t 1 a must be the following: 

(a) For y 0, there are two unbounded solitons of types \j2,3f\ and 1]. The three unbounded 
regions at y ^ are labeled by Ej 1 j 2 ,Ej 1 j e , and Ej e _ lt j e , from left to right. 

(b) For y <C 0, there are / — 2 unbounded solitons of types [j r , jr+2] , where 1 < r < I — 2. The I — 1 
unbounded regions at y <C are labeled by Ej lt j 2 , Ej 2 j 3 , . . . ,Ej s j a+1 , . . . ,Ej l _ 1 j l from left to 
right. 

When we add the index m, going from a soliton graph G for (Gr2/)>o to another one G' for 
(Gr 2i £+i)>o, the unbounded regions at y <§; are labeled by 




Figure 36. Adding the index m to the previous graph, we get a new bounded region 
labeled Ej s j s+ i. 

We now analyze the [j s _i,j s +i] soliton in G' . It was unbounded in G but becomes bounded in G", 
so it must have a vertex v at the bottom of it. Without loss of generality v is trivalent, and because 
GcG', the other two solitons emanating from v must go down from v. By the resonance property, these 
two solitons must be of types [j s -i, %] (on the left) and [x, j s +i] (on the right), where j s _i < x < j s +i- 
The only possibility is x — j s or x — m. We cannot have x = j s because then we would have the 
region Ej s l j s next to a line-soliton of type \j s -x,j s ]. Therefore the soliton [j s -i, js+i] splits into two 
solitons of types [j s _i,m] and [m,j s+ i\. Similarly, in G' the [j s , js+2] soliton splits downward into two 
solitons [j s ,m] (on the left) and [m,j s+2 ] (on the right). It is then easy to see that the line-solitons of 
types [to, js+i] and [j s ,to] intersect at a trivalent vertex (which has two edges emanating upwards from 
it, so is colored black), and the third soliton emanating from that vertex is the soliton [j s ,j s +i]- This 
produces the configuration shown in Figure [36] □ 
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This lemma immediately implies the following. 



Corollary 12.2. Let I be a flag of exponentials for a contour plot Ct (u a) 
mines the soliton graph of the contour plot. 



Then X completely deter- 



Remark 12.6. Note that in the procedure of Proposition 12.1 if one starts with indices It = {j\ < ji < 



< is < 3s+i < ••• < ji—i < ji — n} and adds a new index m such that j s < m < j s +i, then the 



effect on the soliton graph is to make bounded a previously unbounded region E 



3b,3s+i j 



and E, 



new unbounded regions Ej 

Ie+2, ■ ■ ■ ,I n will be labeled E a ^, where both a and b must belong to the set {1, 2, 
belong to {to, to + 1, to + 2, ... , n}. 



and to add two 

It follows that every subsequent region we add by using 

, to}, or both must 



We now associate a triangulation T(X) with every flag of exponentials X. We will show that two flags 
of exponentials X and X' produce the same soliton graph if and only if T(X) and T(I') coincide. 

Corollary 12.3. Consider a soliton graph G for (Gr2„)>o with a flag of exponentials X, and let S = 
{E rs } be the collection of all dominant exponentials labeling the regions of the soliton graph. Then the 
indices appearing in S label the edges and diagonals of a triangulation of an n-gon. We use T(G) = T{X) 
to denote that triangulation. 

Proof. Consider an n-gon with vertices labeled 1,2, ... ,n in order. Any triangulation will consist 
of n boundary edges labeled {1, 2}, {2, 3}, . . . , {n — 1, n}, {1, n}, plus n — 3 diagonals. Now note that 
if we follow the procedure of Proposition |12.1[ then the soliton graph from I3 will have three regions, 
and as we go to 1$, . . . , /„, each step will add one bounded and one unbounded region to the graph. 
Therefore the resulting soliton graph will have 2n — 3 regions, of which n are unbounded (and labeled 
Ei2, . . . , E n -i tn , Ei t n). Finally, note that the procedure of Proposition 12.1 guarantees that we will 
never have two regions E ac and E^d of the resulting soliton graph such that a < b < c < d. Therefore if 
we map each region E a b to the corresponding diagonal between vertices a and b of the n-gon, then the 
resulting collection of diagonals will be a maximal collection of pairwise noncrossing diagonals. □ 



Corollary 12.4. Given a flag of exponentials X, the soliton graph constructed in Proposition 12.1 is 
precisely the one that Algorithm 12.1 associates to the triangulation T(X). In particular, the soliton 



graph depends only on the triangulation T{I). 



Proof. We prove Corollary |12.4| by induction. Suppose it is true for flags of exponentials of length 
I. Then if we follow the steps of Proposition |12.1[ the soliton graph G associated to Ig is precisely 
the one that Algorithm 12.1 associates to T = T{I^ C • •• C It). At this step, T is a triangulation 
of an n-gon whose vertices are labeled 31,32, ■ ■ ■ ,j s ,j s +i, ■ ■ ■ , ji = n in counterclockwise order. Let us 
imagine that these vertices are placed around a circle. Now add a new vertex to on the circle, between 
j s and js+i, together with a triangle t that connects vertices j s , j s +i, and to. This produces a new 
triangulation T' . If we apply Algorithm 12.1 to T' , then the effect on the soliton graph is to add a 
new black vertex v in the center of t, together with three line-solitons to vertices j s ,to, j s+ i labeled 
[to, j s +i], [j s 1 is+i]) [js j m] . These new line-solitons separate three regions corresponding to the diagonals 
of t: the regions are labeled Ej 3t j B+1 (which was already a region, but used to be unbounded), Ej 3t7n , 



and E, 



Let G' denote this soliton graph. 



See Figure 37 and compare it to Figure 36 



It's clear that G' is constructed from G in exactly 
the same way that Proposition |12.1| builds a new soliton graph for I\ C ■ • • C Ie+i from that for 
h C •• • C h, and that T = T(G") = T(h C • • • C h+i). □ 



Finally we are ready to prove Theorem |12.1 



Proof. Recall that to any contour plot we can associate a flag of exponentials X. Proposition |12.1| 
shows precisely how X allows us to reconstruct the soliton graph from X. Corollary |12.3| associates 
a triangulation T(X) to each X, and Corollary 12.4 explains that the soliton graph associated to X is 
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FIGURE 37. Adding a new vertex m to the triangulation together with one new trian- 
gle; and the effect on the soliton graph produced from Algorithm 1 1 2. 1| 



precisely the one which Algorithm 12.1 associates to T(X). Therefore every soliton graph for (Gr2, n ] 



>o 



must be one of those produced by Algorithm |12.1| 



Finally, Lemma 12.2 guarantees that every sequence of subsets I3 C I 4, C • • • C I n = [n] is a flag 



for some soliton graph for (Gr 2 , rl )>o, so every graph produced by Algorithm 12.1 can be realized as a 
soliton graph for (Gr2, n )>o- This completes the proof. □ 
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